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A quantitative evaluation of dilepton sources in heavy-ion reactions is performed taking into 
account both thermal and non-thermal production mechanisms. The hadronic thermal emission 
rate is based on an electromagnetic current-correlation function with a low-mass region (LMR, 
M < 1 GeV) dominated by vector mesons {p, ui, (j)) and an intermediate-mass region (IMR, 1 GeV < 
M < 3 GeV) characterized by (the onset of) a multi-meson continuum. A convolution of the emission 
rates over a thermal fireball expansion results in good agreement with experiment in the low-mass 
spectra, confirming the predicted broadening of the p meson in hadronic matter in connection with 
the prevalence of baryon-induced medium effects. The absolute magnitude of the LMR excess is 
mostly controlled by the fireball lifetime, which in turn leads to a consistent explanation of the 
dilepton excess in the IMR in terms of thermal radiation. The analysis of experimental transverse- 
momentum [qr) spectra reveals discrepancies with thermal emission for qt 1 GeV in noncentral 
In-In collisions, which we address by extending our calculations by: (i) a refined treatment of p 
decays at thermal freezeout, (ii) primordially produced p's subject to energy-loss, (iii) Drell-Yan 
annihilation, and (iv) thermal radiation from t-channel meson exchange processes. We investigate 
the sensitivity of dilepton spectra to the critical temperature and hadro-chemical freezeout of the 
fireball. The p broadening in the LMR turns out to be robust, while in the IMR Quark-Gluon 
Plasma radiation is moderate unless the critical temperature is rather low. 

PACS numbers: 
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I. INTRODUCTION 

Dilepton invariant-mass spectra provide the unique op- 
portunity to directly probe the electromagnetic (e.m.) 
spectral function of the hot and dense medium created 
in energetic collisions of heavy nuclei [H, 0, Q. In the 
low-mass region (LMR, M < 1 GeV), this enables the 
study of modifications of the light vector mesons (V=p, 
u), (j)) caused by their interactions with the surround- 
ing matter particles and/or changes in the underlying 
condensate structure. In the intermediate-mass region 
(IMR, 1 GeV < A/ < 3 GeV), electromagnetic emission 
is expected to become continuum-like with rather well- 
defined strength rendering it a suitable tool to infer tem- 
peratures of the excited system (well) before interactions 
cease ("thermal freezeout"). 

Dilepton measurements in heavy-ion collisions at the 
CERN-SPS (at center-of-mass energies of y^=17.3 and 
8.8 AGeV) have proven the presence of substantial ex- 
cess radiation beyond e.m. final-state decays of produced 
hadrons in both the LMR and IMR [1, 0,11. This, 
in particular, corroborated the presence of interacting 
matter over a duration of ~10-15 fm/c in central Pb- 
Au collisions. Moreover, the spectral shape of the ex- 
cess in the LMR could only be accounted for if major 
medium modifications of the vacuum p-meson line shape 
were incorporated [H, 0, Q . Both the implementation of a 
dropping mass or a strong broadening of its width as fol- 
lowing from hadronic many-body approaches, were com- 
patible with the large enhancement observed at masses 
below the free p mass. In addition, with the same un- 
derlying fireball model (lifetime and temperature evolu- 
tion), reasonable agreement with the enhancement mea- 



sured in the IMR 0, [1] has been established [1] (as well 
as with direct photon spectra [l^, [HI)- The excess for 
M > 1.5 GeV was largely attributed to four-pion type 
annihilations in the hadronic phase, with a strength de- 
termined by the vacuum e.m. spectral function, while the 
calculated Quark-Gluon Plasma (QGP) contribution to 
the thermal yield amounted to ^-^30% (^^ 10-20% in the 
LMR [12). 

A recent, substantial, improvement in statistics and 
mass resolution in low-mass dimuon spectra in In-In col- 
lisions [Tsj shows good agreement with predictions for 
thermal radiation with in-medium p-meson broadening as 
following from hadronic many-body calculations (l^ . [3] . 
The shape of the excess radiation is well described from 
threshold (M = 2m^) to M ~ 0.9 GeV, while the abso- 
lute yield is overpredicted by about 30%. The latter can 
be accommodated by a minor adjustment in the thermal 
fireball evolution (amounting to a 30% lower fireball life- 
time), and after inclusion of in- medium lo and (j) decays 
(whose contribution is mostly localized around their free 
mass), as well as four-pion type annihilation (which sets 
in at masses M > 0.9 GeV), a quantitative description 
of the NA60 invariant-mass spectra from threshold to 
^^1.4 GeV in central In-In collisions emerges [l^ (see also 
Refs. [H, While the inclusive mass spectra are well 
described, the comparison of the calculations with newly 
released transverse pair-momentum (qt) spectra [ll re- 
veals some discrepancies at qx > 1 GeV in semicentral 
In-In. 

In the present article we reiterate the main points of 
our previous study of dilepton invariant-mass spectra at 
the SPS; and extend the analysis to qr spectra. In par- 
ticular, we conduct a detailed analysis of sources at high 
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qt, in terms of (i) an improved treatment of p decays 
at thermal freezeout (which are subject to maximal blue 
shift due to transverse flow), (ii) a component of pri- 
mordial (hard-produced) p mesons subject to energy loss 
when traversing the medium (using high-px pion spec- 
tra as a guideline), (iii) Drell-Yan (DY) annihilation in 
primordial N-N collisions which we extrapolate to small 
mass by imposing constraints from real photon produc- 
tion, and (iv) meson ^-channel exchange contributions 
to the thermal production rate which are not included 
in the many-body vector-meson spectral functions (but 
which have been found to be a significant source at high 



qx in real photon production [ll|) . Another interest- 
ing issue which has received little attention thus far is 
how uncertainties in the critical temperature and hadro- 
chcmical evolution of the fireball affect dilepton spectra. 
For hadro-chcmical freezeout we investigate the sensitiv- 
ity to temperatures in the range Tch — 160-175 MeV, rep- 
resentative for top SPS energy according to recent ther- 
mal model analyses [2l|, . In connection with updates 
of lattice-QCD results indicating a critical temperature 
up to Tc ~ 190-200 MeV [1^, this, in particular, opens 
the possibility of a chemically equilibrated hot and dense 
hadronic phase for, say, T = 160-190 MeV, which we 
also consider. We furthermore conduct a quantitative 
study of effective slope parameters of the NA60 qr spec- 
tra, where the investigations of the fireball chemistry are 
supplemented with variations of the transverse flow ve- 
locity. Finally, wc revisit the consequences of our fireball 
refinements on our previous evaluations (23 | of diclcc- 
tron spectra as measured by CERES/NA45 in semicen- 
tral Pb-Au collisions [1] , as well as recent data in central 

Pb-Au m. 

Our article is organized as follows: in Sec. [II] we recall 
the main ingredients of our approach to calculate thermal 
emission rates, based on the e.m. spectral function in the 
vacuum fSec. lII A]) , followed by discussing medium effects 
on hadronic emission at low mass (due to in-medium p, 
Lo and (j) spectral functions; Sec. Ill Bp and at intermedi- 
ate mass (due to finite-T chiral mixing; Sec. Ill C"l) : a new 
element not included in previous spectral-function cal- 
culations are ^-channel meson-exchange reactions which 
therefore are elaborated in more detail in Sec. Ill Dl while 
we will be brief on partonic emission from a Quark-Gluon 
Plasma (Sec. IIIEp . In Sec. IIIII we evaluate nonthermal 
dilepton sources: p-meson decays at thermal freezeout 
(whose decay kinematics differ from thermal radiation; 
Sec. nil Ap . an estimate of primordial p's at large qr 
(which do not thcrmalize; Sec. IIIIBj) . and Drell-Yan an- 



^ Note that contributions {ii)-{iv) have little bearing on the inclu- 
sive invariant-mass spectra which are predominantly populated 
by low-momentum sources with ^ 1 GeV. Therefore, the in- 
clusion of these contributions does not upset our earlier descrip- 
tion of inclusive M spectra in terms of thermal radiation and 
freezeout p's. Initial results of these studies have been reported 
in Refs. [lllli. 



nihilation (with an extrapolation to low mass; Scc. lIII C|) . 
In Scc. lIVI we recollect the ingredients to our thermal fire- 
ball model for the space-time evolution of the medium in 
heavy- ion collisions, including variations in hadrochem- 
ical freezeout and critical temperature. In Sec. |V] we 
implement all dilepton sources into the fireball to com- 
pute dimuon invariant-mass (Sec. IV Ap and transverse- 
momentum spectra (Sec. IV Bp in comparison to NA60 
data; the effect of hadrochemistry and Tc on NA60 spec- 
tra is worked out in Sec. lVC| followed by a slope analysis 
of qr spectra (Sec. IV Pp in the context of which we also 
investigate different radial fiow scenarios; the improve- 
ments in the fireball and dilepton source description as 
deduced from the NA60 data are confronted with previ- 
ous and new CERES/NA45 data in Sec. EH Sec. EJ 
contains a summary and conclusions. 



II. THERMAL DILEPTON RADIATION 



A. Emission Rate and Electromagnetic Spectral 
Function 



In thermal equilibrium the rate of dilepton emission 
per four-volume and four-momentum can be related to 
the hadronic e.m. spectral function as [1^ 



(1) 



X /^(go;T) , 



which in this article we will refer to as thermal dilep- 
tons (or thermal radiation). The retarded e.m. current- 
current correlator is given by 

ne^^(g) = i / d^xc'«'-'e(x°) ([J^:„(a:), J:^(0)]) , (2) 



where a=e^/(47r)=l/137 denotes the fine structure con- 



stant. 



q the dilepton invariant mass squared 



with energy go a-nd three-momentum q, and f^{qQ;T) 
the thermal Bose distribution function (T: temperature, 
Pb- baryon chemical potential). The final-state lepton 
phase space factor. 



L{M) = 1 



2m[ 



1 - 



4mf 
M2" 



(3) 



depends on the lepton mass, m;=0. 511(105. 6) MeV for 
electrons (muons; I = e,^), but quickly approaches one 
above threshold, M = 2mi {e.g., L{M = 0.3 GeV) = 0.89 
for dimuons). 

In the vacuum, the e.m. spectral function, Imneni(-^), 
is well known from e^e~ annihilation into hadrons. It 
is characterized by the light vector resonances p(770), 
a;(782) and (/)(1020) at low mass (vector dominance model 
(VDM)) and a pcrturbative quark-antiquark continuum 
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at higher mass, 



E 



V—p,Lij,(i> 

12-n 



av 



ImDv. M < M, 



dual 



(1 



lE(e.)^ M>M, 



dual 



(4) 

{Nc = 3: number of colors, a^: strong couphng con- 
stant, Ci'. the electric quark charge in units of the elec- 
tron charge, and i is running over up, down and strange 
quark flavors for M < 3 GeV). Mduai — 1-5 GeV signifies 
a "duality" scale above which the total cross section for 
e+e^ — > hadrons (and thus the strength of Imllom) es- 
sentially behaves perturbatively with little impact from 
subsequent hadronization. Since the hadronic final state 
in e"'"e~ annihilation approximately resembles a thermal 
medium (except for strangeness), time- reversal invari- 
ance of strong and electromagnetic interactions implies 
that, to leading order in temperature, the equilibrium 
dilepton emission rate of hadronic matter is determined 
by the free e.m. correlator, n™^, cf. also Refs. (ol. [27|. 

In the hadronic basis, the e.m. spectral function is 
dominated by the isovector (p) channel while the isoscalar 
channels are suppressed. The SU(3)-flavor quark model, 
e.g., predicts a weighting of 9:1:2 for the e.m. cou- 
plings of p, Lo and (j), respectively. These values are 
roughly in line with the electromagnetic decay widths 
Tv^ee = ^mv/gy = 7.0, 0.60 and 1.27 keV for 
p, uj and (f), respectively (the empirical values of the 
width should be corrected for phase space by the ratio 
m^/m0=O.77). 

The vector-isovector channel furthermore provides the 
most direct link to chiral symmetry: under chiral ro- 
tations, the p channel transforms into the axialvector- 
isovector (ai) one. Thus, the ai is commonly identified 
as the chiral partner of the p.^ The pertinent spectral 
functions in free space have been accurately measured 
in hadronic t decays H^, H^l, cf. Fig. [U and evaluated 
in terms of chiral order parameters (/jr, four-quark con- 
densates) using Weinberg sum rules [3l|, [13] ■ To leading 
order in temperature {i.e., for a hot pion gas), the mu- 
tual in-medium modifications of vector and axialvector 
channels can be inferred from chiral symmetry alone, as 
will be discussed in Sec. Ill CI below. 

The main objective for the remainder of this Section 
is a realistic evaluation of the in-medium e.m. spectral 
function. In the LMR, we assume that vector domi- 
nance remains valid in the medium. Our calculations 
of the F-meson spectral functions. Ay = —2lmDv, 
utilize effective hadronic Lagrangians, where the inter- 
action vertices and coupling constants are constrained 



^ An alternative has recently been suggested in Ref. [2^ 1 in terms 
of the "Vector Manifestation" (VM) of chiral symmetry, where 
the chiral partner of the (longitudinal) p is identified with the 
pion. 
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FIG. 1: (Color online) Free isovector-vector (upper panel) 
and -axialvector (lower panel) spectral functions as measured 
in hadronic r-decays [2y|, compared to fits with three- and 
four-pion contributions (the two-pion piece follows from the 
previously calculated p propagator [3^). 



by gauge and chiral symmetry, as well as empirical de- 
cay and scattering data. Medium effects are calculated 
within hadronic many-body theory without introducing 
explicit medium dependencies of the bare parameters in 
the Lagrangian (which require information beyond the 
effective hadronic theory). A careful comparison of our 
predictions for dilepton spectra in heavy-ion collisions 
with experiment can then serve as a basis for identifying 
medium effects that go beyond the hadronic many-body 
framework. 



B. Hadronic Emission at Low Mass: 
Vector Mesons 



In-Medium 



Since the prevalent role in the LMR of the e.m. cor- 
relator is played by the p meson, the latter has been 
the main focus for studying in-medium effects in thermal 
dilepton production. However, with the current precision 
of the NA60 dilepton data on the 10-20% level, thermal 
emission from tu and (j) decays becomes relevant. Their 
contributions are an inevitable consequence of the forma- 
tion of a thermal medium, and we therefore incorporate 
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them along with pertinent medium modifications. 

Within the VDM the in-medium e.m. correlator in the 
LMR remains directly proportional to the vector-meson 
spectral functions as in Eq. (jlj. The key objective is 
then to calculate the vector-meson selfenergies at finite 
temperature and baryon density (qb), Sy, figuring into 
the propagator as 



Here, rriy^ denotes the bare mass, and the selfenergy 
contributions are classified into three types: (i) medium 
modifications of the pseudoscalar meson cloud, Syp 
{P = TTTT, Stt or KK for V = p, uj or (p, respectively), 
and direct interactions of V with (ii) mesons and (iii) 
baryons from the surrounding heat bath, E^a/ and Sys. 
One should also note that without the implementation of 
a phase transition into a realistic hadronic model (which 
has not been achieved thus far) , many-body calculations 
are to be considered as an extrapolation beyond a cer- 
tain temperature and/or density. We estimate that, un- 
der SPS conditions, our approach should be reliable up 
to temperatures of about T ~ 150 MeV, where the to- 
tal hadron density amounts to about — 2po (with 
a net-baryon density of about O-lgo), while it has in- 
creased to about 5po a-t the expected phase boundary at 
Tc ~ 175 MeV. In the following we briefly recall the con- 
struction of, and constraints on, the underlying effective 
Lagrangians. 



1. p Meson 

The first step in the calculation of a realistic p propa- 
gator is a proper description of its coupling to the pion 
cloud in vacuum (giving rise to its two-pion decay width, 
rp7rir(wp) ~ 150 MeV). A suitable chiral efi^ective La- 
grangian can be obtained by, e.g., introducing the vec- 
tor fields as local gauge fields into the chiral pion La- 
grangians and implementing the photon field via vector 
dominance (3^ . ISq , 

i^-npi =9pPp ■ (tT X 5^7?) 

+ -^A,pt^ . 
9p 

The pertinent two-pion and one-pion-tadpole diagrams 
specify the p's pion cloud in the vacuum with a transverse 
selfenergy (using a proper regularization procedure). The 
coupling constant, gp, a cutoff parameter, (to ren- 
der the pion loops finite), and the bare p-mass, mp°\ 
can be readily adjusted to reproduce P-wave tt-tt scat- 
tering phase shifts and the pion e.m. formfactor in free 



space [33[ . This model also predicts the two-pion part of 
the r-decay data well, cf. Fig. [T] 

Based on the above model in the vacuum, we adopt 
the in-medium p propagator from the hadronic many- 
body approach developed in Refs. [H, [H, [s^, to which 
we refer the reader for more details. Modifications of 
the pion cloud are induced by an in-medium tt propa- 
gator which is dressed with standard NN~-^ and AN~^ 
excitations in nuclear matter, supplemented with appro- 
priate vertex corrections to maintain a conserved vec- 
tor current. While the calculations of nuclear pion-cloud 
modifications are often restricted to the case of vanishing 
p three- momentum ("back-to-back" kinematics), we here 
employ the extension to finite three-momentum worked 
out in Ref. [ssj . This is particularly important in that it 
enables (i) to constrain the corresponding medium effects 
by nuclear photo-absorption data (sj and nN pN 
scattering, and (ii) a meaningful application to dilepton 
transverse momentum spectra. At finite temperature, 
where a (large) fraction of nucleons is thermally excited 
into baryons. only the nucleon density should be em- 
ployed. To account for (i) ttBN~^ excitations into higher 
resonances and (ii) nB2B^^ excitations (Bi ^ N), an 
"effective" nucleon density g'j^ = g^ + O.Sps {B ^ N) 
has been used in the pion propagators [24| (we estimate 
the uncertainty implied by this prescription to be about 
±10%). In addition, effects of antibaryons are added by 
multiplying with a factor (1 -I- p/p) where p/p de- 
notes the measured antiproton-to-proton ratio (this is a 
small correction at SPS energies but becomes essential at 
RHIC nil ) . The modifications of the pion cloud induced 
by a pion gas have been found to be rather small in vari- 
ous analyses [H SO, [HI, S El . A large part of the pion- 
gas effect is already provided by the Bose enhancement 
factor, [1 -I- /'^(wi) -I- f'^{LU2)] [41[, while low-energy tt-tt 
interactions are suppressed by chiral symmetry. The for- 
mer is included in the present calculations, but the latter 
are neglected (being much weaker than baryon-induced 
effects). 

Direct interactions of the (bare) p with baryons and 
mesons have been approximated by resonance interac- 
tions (pion-cxchangc interactions, including spacelikc pi- 
ous, are implicit in the nuclear modifications of the pion- 
cloud part). For p + N — > B {i.e., pBN~^ excita- 
tions), this includes both 5- wave scattering into negative- 
parity baryons and P-wave scattering into positive- 
parity baryons, with B = iV(1520), A(1620), A(1700) 
and N, A(1232), iV(1440), 7V(1720), A(1905), 7V(2000), 
respectively. An estimate of the involved couplings and 
cutoff parameters can be obtained by reproducing the 
empirical B pN, jN decay widths, but more reliable 
constraints follow from comprehensive scattering data. 
For the present model, in addition to hadronic decay 
branchings, cold nuclear matter effects induced in both 
the pion cloud and via resonant p-N interactions have 
been checked against total photoabsorption cross sections 
on the nucleon and nuclei [s^l, as well as nN pN 
scattering data (T], H^] . Whereas 7iV and pN scattering 
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relates to the low-density limit of the p-spectral func- 
tion, J A data constrain lmDp{M 0,q) close to nu- 
clear saturation density, go = 0.16 fm^''. For the finite- 
temperature case, excitations on hyperons of the type 
pY2Yf^ are incorporated for Yi = A(1115), I](1192) and 
Y2 resonances with quantum numbers equivalent to the 
nonstrange sector and appreciable empirical decay widths 
into p's and/or photons. pB2B^^ excitations on ther- 
mally excited nonstrange baryons {e.g., A(1930)A~-'^) are 
included but turned out to be small. 

In the meson gas, direct p-M interactions on ther- 
mal mesons {M = tt, K, p) have been scrutinized in 
Ref. psf . again based on s-channel resonance dominance, 
pM — > R, and with pertinent Lagrangians satisfying ba- 
sic requirements of vector-current conservation and chiral 
symmetry. As in the case of excited baryons, the cou- 
plings and cutoffs have been estimated from combined 
hadronic and radiative decay widths. The most impor- 
tant contributions to the p propagator arise from axi- 
alvector mesons, R — ai(1260) (chiral partner of the /?), 
/ii(1170) and ivri(1270), as weh as the w(782). 

Another valuable consistency check for the in-medium 
properties of hadrons arc QCD sum rules. The latter 
relate a dispersion integral over a spectral function to 
an operator product expansion (OPE) for the correla- 
tion function at spacelike four-momenta, f/^q'' < 0. The 
nonperturbative coefficients of the OPE involve quark 
and gluon condensates whose in-medium modifications 
reflect upon (an energy integral over) the spectral func- 
tion. Generically, a reduction of the condensates leads 
to a softening of the spectral function, i.e., a shift of 
strength to lower mass. However, as has been elaborated 
in Ref. [1^1 for the case of cold nuclear matter, the ad- 
ditional low-mass strength in the in-medium p spectral 
function can be due to both a reduced mass or a broaden- 
ing of the width, or a combination thereof. The pertinent 
"allowed band" for width and mass values at normal nu- 
clear matter density has been quantified in Ref. (46j . For 
the p spectral function employed here [l^l, one finds a 
mass and width of (mp,Vp)\f,^ ~ (820,450) MeV at nu- 
clear saturation density and vanishing three-momentum. 
These values turn out to be consistent with the band de- 
rived in Ref. [1^ (cf. middle panel of Fig. 2 in there), 
but more extensive (i.e., different densities and temper- 
atures) and quantitative comparisons are certainly desir- 
able (e.5., the low-mass strength (shoulder) in the p spec- 
tral function which is induced by "Pisobar" (ttAA^"^) 
and "Rhosobar" {pBN~^) excitations, cannot be rep- 
resented in Breit-Wigner form as underlying Ref. [3). 
However, at this point, there is no indication for changes 
in the bare parameters in the Lagrangian, i.e., at nuclear 
saturation density hadronic many-body effects saturate 
the reduction of quark and gluon condensates as required 
by QCD sum rules. 

The p-spectral functions under fireball conditions rep- 
resenting the evolution of the hadronic phase at full SPS 
energy {y/s = 17.3 AGcV) arc displayed in the upper 
panel of Fig. [21 The resonance strongly broadens with lit- 
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FIG. 2: (Color online) In-medium spectral functions of p 
(upper panel), lo (middle panel) and </> (lower panel) under 
conditions representative for the time evolution of heavy- 
ion collisions at full SPS energy: (pjv , p-n , Pk , T) [MeV] = 
(232,0,0,175), (331,39,62,150), (449,79,136,120) and vac- 
uum for the dotted, short-dashed, long-dashed and solid lines, 
respectively. 



tie, if any, mass shift, essentially leading to its "melting" 
close to the expected phase boundary. These features are 
a rather generic consequence of a multilevel excitation 
spectrum in which real parts occur with altering signs 
while imaginary parts are negative definite and strictly 
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add up. Also note the development of a pronounced low- 
energy strength around M = 0.4 GeV which leads to a 
very large enhancement over the vacuum spectral func- 
tion (and will be augmented by thermal Bose factors in 
the dilcpton rate). At lower temperatures, the presence 
of meson-chemical potentials (for tt, K , rf) implies higher 
meson and baryon densities than in chemical equilibrium 
at a given temperature, entailing stronger medium ef- 
fects. The build-up of the chemical potentials (necessary 
to preserve the experimentally observed hadron ratios 
throughout the hadronic cooling) therefore facilitates the 
observation of in-medium effects in the dilepton spectra. 
A strong broadening of the p in hot and dense matter 
has also been reported in other calculations at finite T 
and |43,|ii. 

The in-medium p-spectral function has been rather ex- 
tensively employed to calculate dilepton spectra in heavy- 
ion collisions. At full SPS energy it has been imple- 
mented into fireball models [U, [sHi, transport simula- 
tions [4^ and hydrodynamic evolutions [50], with fair 
success (and consistency) in describing CERES/NA45 
data in central and semicentral Pb-Au collisions [1]. 
When averaged over a typical space-time evolution, the 
in-medium p width at SPS amounts to r™^'^ ~ 350- 
400 MeV, almost three times its vacuum value, implying 
T piTc) — nip when extrapolated to the expected QCD- 
phasc boundary. Despite an experimental pion-to-baryon 
ratio of ^ 5 : 1, medium effects induced by baryons have 
been identified as essential at SPS energies early on [3^. 
This has led to the prediction of an even larger excess 
at lower beam energies, which was tentatively confirmed 
in a 40 AGeV Ph+Au run of CERES/NA45 H (albeit 
with rather large errors). The lightlike limit of the p- 
spectral function has been applied to calculate thermal 
photon spectra supplemented with meson-exchange 
reactions which are not included in Im Dp (these become 
significant at momenta above 1 GeV and will be discussed 
in Sec. Ill Dl below). After convolution over the same fire- 
ball model as used for dileptons, the direct-photon excess 
observed by WA98 in central Pb-Pb reactions at SPS p^ . 
as well as by PHENIX in central Au-Au collisions at 
RHIC [5l| , can be reasonably well accounted for (with ad- 
ditional primordial and QGP photon sources, the latter 
being subdominant at SPS). The in-medium p-spectral 
function following from hadronic many-body theory 
has thus passed a rather wide range of both theoretical 
and phenomenological tests, in both cold nuclear as well 
as hot hadronic matter. 



2. uj Meson 

Medium effects on the u) have received less attention so 
far, especially in the context of ultrarelativistic heavy- ion 
collisions (URHICs) where its dilepton decays are domi- 
nated by free decays after thermal freezeout (commonly 
referred to as a "cocktail" contribution). We here adopt 
the same approach as for the p, following Ref. [s^. In 



the vacuum, the hadronic lo width, T^^^t,, is composed 
to approximately equal parts of a direct cj ^ Stt and a 
Lj ^ pn coupling, represented by anomalous terms in the 
interaction Lagrangian, 

CujpTT = G^pT, epvar d^p" ■ 9^7? , (7) 
C^z. = G^3. ^^...r c^'' (a'' 7? X a-7f) • a" 7? . (8) 

With a hadronic formfactor cutoff-parameter, A^p^: = 
1 GeV, and using the VDM (4^, one finds for the ra- 
diative decay width V^^^-^-^ = 0.72 MeV which approx- 
imately agrees with the updated experimental value of 
0.76 ±0.03 MeV 

In cold nuclear matter, the renormalization of the n-p 
cloud has been found to be rather sensitive to the medium 
effects on the p (ssl . IZsl . [ssl . [53 |. which is due to the 
fact that the oj Trp{—> tttttt) decay proceeds via the 
low-mass tail of the p-spcctral function below its nomi- 
nal mass of 770 MeV. In Ref. ^ 80% of the estimated 
in-medium broadening at nuclear saturation density is 
attributed to the in-medium p spectral function (taken 
from Ref. [s^l), resulting in Tujpigo) — 60 MeV. Like- 
wise, a large contribution to second order in the nuclear 
density, 0{g%), to FupiQN) arises due to the simultane- 
ous dressing of both tt and p [11] . Following Ref. , we 
restrict ourselves to the impact of the in-mcdium p by 
folding the phase space in the w — s- 7rp-decay width over 
the spectral function evaluated in the previous section, 

J\fmax 

2Glp^ f MdM 


X [1 + r{u;^) + fP{Ep)] F^p^{q,^f , 

where M^^^=y/s-m^, ul = ml + q^^, = + q^^, 
and ^cm is the three-momentum of tt and p in the rest 
frame of the lu. In addition we account for cold- nuclear 
matter effects due to the two most important direct uj-N 
excitations, ujN{1520)N~'^ and ivN {1650) N'^ , as iden- 
tified in Ref. (ssj in a coupled-channel analysis of empir- 
ical vector-meson nucleon scattering amplitudes. Recent 
measurements of 7r*'7-invariant mass spectra in photon- 
induced production off hydro gen and Nb targets have 
revealed a low-mass shoulder [Sa] that has been inter- 
preted as reduced lo mass in nuclear matter. However, 
this interpretation is not free of controversy and may also 
be explained by a large increase of the uj width to about 
90 MeV (53. The latter is compatible with the uj prop- 
agator employed in the present paper. 

Finite-temperature effects on the lo propagator in- 
clude Bose-enhancement factors in the n-p cloud, the in- 
medium p propagator in Eq. a schematic estimate of 
the inelastic ttuj — > tttt width based on Ref. [sst . and a 
full calculation of the thermal w -I- tt 6i(1235) loop [38| . 

The predicted average lo width in the hadronic phase of 
URHICs is r^<='i~100 MeV [H, cf. middle panel of Fig.d 
quite consistent with the results of Ref. [53. This has 
the interesting consequence [11] that if dilepton spectra 
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in HICs show a structure of '--'100 MeV width in the p-joj 
mass region, it would be an unambiguous signature of the 
in-medium w spectral function since it is narrower than 
the lower limit given by a free p contribution, F™""^ < 
~ 150 MeV. 

3. (j) Meson 

For the 0, collision rates in a meson gas have been 
estimated to amount to a broadening by ^ 20 MeV at 
T = 150 MeV [gOI- The dressing of the kaon cloud is pre- 
sumably the main effect for (j) modifications in nuclear 
matter, increasing its width by ~ 25 MeV at satura- 
tion density [6l| . Recent data on absorption in nuclear 
photoproduction suggest substantially larger values [6^ . 
Since a quantitative, empirically constrained treatment 
of the (j) in hot and dense matter is not available at 
present, we will consider the <j) width as a parameter. Fol- 
lowing Ref. [iBl we augment the microscopic calculations 
of Ref. [11] by a factor of 4 to roughly agree with the 
phenomenological values extracted from nuclear photo- 
production |62l |. The corresponding spectral functions 
under SPS conditions arc displayed in the lower panel 
of Fig. [21 The average width over the fireball evolution 
amounts to F'°°^ ~ 80 MeV. 



C. Hadronic Emission at Intermediate Mass: 
Multi-Pion Annihilation 

Above the </> mass, the hadronic structure of the (vac- 
uum) e.m. correlator becomes more involved being char- 
acterized by overlapping broad resonances which com- 
bine into continuum-like multi-meson states (most no- 
tably isovector four-pion states), recall Fig. [T] To es- 
timate medium effects in this regime, we take recourse 
to model-independent predictions following from a low- 
temperature expansion evaluated using chiral-rcduction 
techniques for the corresponding matrix elements involv- 
ing thermal pions. This method has been first worked 
out in Ref. |63| in the chiral limit (m^ = 0) and leads to 
the chiral mixing formula for the isovector part of vector 
and axialvector correlators, 

ny(g) = (i-£)nr(g)+enr(g), 

n^(<?) = enr(g) + (i-e)nr(9). ^ ^ 

The mixing coefficient e is determined by pion tadpole 
diagrams via a loop integral 

{ijj^: on-shell pion energy, /^r = 93 MeV: pion decay con- 
stant). For £ — > 1/2, V and A correlators degenerate sig- 
naling chiral-symmctry restoration (at T^'^cilGO MeV). 
Note that the admixture of the axialvector part in the 



M = 1-1.5 GeV region (which corresponds to the oi res- 
onance) fills in the "dip" structure of the vector corre- 
lator (recall Fig. [T]); for e ^ 1/2 both spectral func- 
tions merge into the perturbative-QCD continuum for 
M > 1 GeV which has been interpreted [l| as a lower- 
ing of the quark-hadron duality scale, Mduai, cf. Eq. jl]). 
The simple form of Eq. (fTO|) only holds in the soft-pion 
limit, i.e., when neglecting the pion four-momenta in 
Ily^^. More elaborate treatments [H, [6^ will broaden 
and somewhat reduce the enhancement in the M ~ 1- 
1.5 GeV region [b^. We here implement the mixing ef- 
fect on the (isovector) four-pion part utihzing Eq. ([TU]) 
with a mixing parameter e = ■^£{T, Pt^) /e{Tc,0) with fi- 
nite pion mass (mir = 139.6 MeV), critical temperature, 
Tc, and pion chemical potentials (p-^ > 0) implemented 
as an overall fugacity factor, = e'*"/^, in Boltzmann 
approximation (pt^ varies according to the thermal fire- 
ball evolution discussed in Sec. IIVI below). The two-pion 
piece, as well as the three-pion piece corresponding to 
ai decay, oi n + p, have been removed as they are 
included via the p-spcctral function as discussed above. 
A detailed analysis of the derivation of Eq. (jlOp in the 
presence of a finite pion-chcmical potential leads to the 
following expression for the mixing effect on the vector- 
isovector current correlation function: 

uyiq) ={1 - e)ztui% + |4nr3. + 1(4 + 4Wa% ■ 

(12) 

D. Hadronic Emission at High Momentum: Meson 
£-ChanneI Exchange 

As discussed in the Introduction, a current discrepancy 
in the theoretical description of the NA60 data concerns 
the dilcpton yield at high transverse pair-momentum, 
qt > I GeV, for masses M < 1 GeV. It is therefore 
important to scrutinize the vector-meson spectral func- 
tion calculations for high momentum sources. The latter 
are typically associated with t-channel meson exchanges, 
rather than s-channcl resonances whose decay products 
are suppressed at momenta far off the resonance. In the 
spectral-function approach, t-channel exchange processes 
are encoded in medium modifications of the pseudoscalar 
meson cloud. As elaborated in Sec. lIIBTl the pion-cloud 
modifications of the p explicitly include interactions with 
baryons, but not with thermal mesons. Consequently, 
t-channel meson exchanges are present for pB interac- 
tions, but not for pn. E.g., the experimental nN — > pN 
cross section is properly reproduced at large y/s where t- 
channel w exchange dominates. Indeed, it has been found 
in Ref. [TU in the context of thermal photon production 
that t-channel exchanges ulp + ^T—^^T + "f outshine the 
contributions from the p-spectral function at momenta 
above ~1 GeV for T = 200 MeV (at smaller tempera- 
tures the importance of t-channel exchanges is reduced) . 
Most notably, uj exchange emerged as the single most 
important hadronic high-^/^ thermal photon source. 
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FIG. 3: The Feynman diagram corresponding to the annihi- 
lation of a p meson via a;-meson exchange in the t channel. 

To estimate the relevance of t-channel exchanges in irp 
scattering for thermal dileptons, we adopt the same sim- 



plified treatment as in Ref. properly extended to a 
finite invariant mass of the virtual photon (lepton pair) : 
rather than implementing the pertinent diagram into the 
p selfenergy, we employ the standard kinetic theory ex- 
pression for the thermal production rate 

(13) 

which allows to calculate the contribution for tt + p — > 
TT + / + + in terms of the corresponding matrix element, 
Mttp^ttU, in Born approximation. The latter substan- 
tially facilitates the task of maintaining gauge invariance. 
For the pujT: vertex, the anomalous-coupling Lagrangian, 
Eq. (O, is used to calculate the matrix element, 



4 (^2 



X {{M^ 



\t- 



;L{M)Q{M - 2mi) 



+ + 2ml + m^pf - 2{M^ + 2ml + + 2s^] - 2t^{M'^ + 2ml + - s) + t^} , 



(14) 



r 



corresponding to the Feynman diagram depicted in 
Fig. [3l M denotes the invariant mass of the lepton pair; 
s = (pi +P2)^{Pi +P2)fi and t ~ k^'kf^ are the usual Man- 
delstam variables with k'^ = {p^ — pi)^ the exchanged 
four-momentum. The p-7 vertex has been written in the 
form [il 



C 



PI 



-Cm'pA^p^" 



(15) 



where C = e/gp — 0.052. Medium effects on the p meson 
are implemented by using the full in-medium propagator, 
Dp{M), for the intermediate p meson which couples to 
the photon (upper p line in Fig.[3l with = g^) and by 
folding Eq. (|14p over the in-medium spectral function for 
the incoming p (lower p line in Fig. [3l with = p^ and 
a weight —2mp/TrlmDp(mp)). To avoid double counting 
with the to s-channel graph already included in the in- 
mcdium p-selfenergy, the integral in Eq. p3p is restricted 
to spacelike w-exchangc momenta, i.e., t < 0. The cou- 
pling constant, ^p^^ = 25.8 GeV~"^, has been fixed in [i^ 
by a simultaneous fit to the hadronic and radiative uj de- 
cays including a hadronic dipole-form factor. 



F{t) 



2A2 



2A2 -i 



(16) 



with A = 1 GeV. To maintain gauge invariance in a sim- 



plified way as in Ref. [ll| , we pull the formfactor outside 
the integral and introduce an average momentum trans- 
fer, t (adapted to the finite mass of the virtual photon), 
according to 



2 ./ \m^. 



{ml - t)2 2 J [ml - t{xW ' 
-1 

t{x) = +ml- 2{M^ + - qpx) . 



(17) 



(18) 



The amplitude, Eq. ([14]), is then multiplied by F^(i). 
Eq. (|13p furthermore accounts for pion fugacity factors, 
Ztt due to finite pion chemical potentials, implemented in 
Boltzmann approximation {e.g., for the incoming tt and p 
this amounts to a total fugacity of z^). In Fig.|4]we sum- 
marize our results for the uj i-channel emission rate, in- 
tegrated over three-momentum, for conditions along our 
default trajectory for In(158 AGeV)+In collisions at SPS 
(including finite fi^^). The upper panel illustrates that the 
implementation of the in-medium p propagator (which 
was not done in Ref. [ll|) leads to a notable reduction 
in the emission rate, which is mostly due to the reduc- 
tion in the intermediate p propagator, \Dp{M)\'^. We 
have verified that medium effects on the uj propagator are 
negligible (in the spacelike regime finite- width effects are 
insignificant as long as F^ ^ m^). In the middle panel. 
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FIG. 4: Three-momentum integrated thermal emission rates 
from Lj t-channel exchange in the vrp — > 7re''"e~ reaction. Up- 
per panel: comparison of the rates when using either the free 
or the in-medium p propagators in both the matrix element 
and the incoming p-mass distribution. Middle panel: compar- 
ison of the in-medium uj f-channel rates with thermal rates 
from the full in-medium p spectral function. Lower panel: 
same as middle panel, but restricting the momentum integra- 
tion to <j > 1 GeV. All curves for r=150 MeV and 120 MeV 
include pertinent pion fugacities at p^ = 39 MeV and 79 MeV, 
respectively. 



the comparison of the in-medium w i-channel rates to 
the ones from the full in-medium p spectral function [2^ ] , 
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FIG. 5: Mass-integrated dilepton rates as a function of three- 
momentum for LO t-channel exchange in -Kp — > 7re^e~ (com- 
puted with in-medium p propagators) and the leading con- 
tribution from the full in-medium-p spectral function. As in 
Fig- HI the results for T=150 MeV include pion fugacity fac- 
tors. 



based on Eq. ([T]), confirms that the former contribution 
is indeed small in magnitude. However, if one applies a 
three-momentum cut of g > 1 GeV (lower panel) , the rel- 
ative magnitude of the f-channel contribution increases, 
as anticipated at the beginning of this section. But even 
in this case the relative strength of the oj ^-channel ex- 
change is rather moderate. The increasing contribution 
at higher momentum can be more directly seen when in- 
tegrating the rate, Eq. (|13p . over invariant mass bins and 
plotting it versus q, cf. Fig. [S] 

As indicated above, in our calculations of dilepton 
spectra in Sec.|Vl we will implement the uj t-channcl emis- 
sion rate, Eq. p3p . by an incoherent addition to the main 
contribution. Eq. ([T|). This neglects interference terms in 
the pertinent selfenergy contribution in the (denomina- 
tor of the) p spectral function, which is justified due to 
the relative smallness of the ui i-channel part. To simu- 
late the presence of other f-channel processes [e.g., tt and 
ai exchange), guided by the photon rate calculations of 
Ref. pdj . we will multiply the w t-channel contribution 
by a factor of 2 in our calculation of dilepton spectra 
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below. 



E. Partonic Emission: Quark-Antiquarlc 
Annihilation 



Emission from the QGP is calculated using the hard- 
thermal-loop improved rate for qq annihilation [67j . in- 
cluding an extrapolation to finite three-momentum. It 
turns out [111 that this rate has the conceptually at- 
tractive feature that it closely coincides with the rate in 
hadronic matter when both are extrapolated to the ex- 
pected phase-transition region. This is suggestive for a 
kind of quark-hadron duality and has the additional 
benefit that the emission from the expanding fireball can 
be anticipated to become rather insensitive to details of 
how the phase transition is implemented {e.g., to the val- 
ues for the critical temperature or "latent heat"). This 
point will be studied explicitly in Sec. IV Cl below. 

Recent (experimental and theoretical) developments 
suggest that the QGP features significant nonperturba- 
tive effects for temperatures up to ^1.5-2Tc, e.g., the sur- 
vival of (hadronic) resonances and/or bound states, even 
in the light-quark sector. Possible consequences for the 
dilepton emission rate have been estimated in Ref. [6^. 
Depending on the width of these states, a maximal en- 
hancement of up to a factor of ~ 2 over the perturbative 
QGP emission rate at intermediate masses {i.e., in the 
mass range of 1.5-2 GeV) is conceivable. We will not 
further pursue this possibility in the present paper. 



III. DILEPTON SOURCES OTHER THAN 
THERMAL RADIATION 

In this section we address dilepton sources other than 
radiation from a thermal source as given by Eq. ([T]), i.e., 
p decays after thermal freezeout, decays of primordially 
produced p mesons which escape the fireball, and pri- 
mordial Drell-Yan annihilation. For simplicity, we generi- 
cally refer to these sources as "non-thermal" , even though 
the freezeout decays, e.g., are represented by thermal 
"blast-wave" spectra commonly used to characterize light 
hadron spectra in heavy-ion collisions. 



A. p Mesons at Thermal Freezeout 

Dilepton decays of long-lived hadrons (most notably 
Dalitz decays of r] and to, as well as exclusive de- 
cays of Lo and 0) can be rather well separated from the 
freezeout of the interacting fireball in heavy-ion collisions, 
which has given rise to the notion of the "hadronic decay 
cocktail" contribution to dilepton spectra, computed us- 
ing the free spectral shape of the decaying mesons. For 
the p meson, due to its short lifetime, such a separation 
is not well defined. Therefore, when adding thermal ra- 
diation to the cocktail, the p is commonly removed from 



the latter and implemented into the thermal yield. In 
previous works [l^, [U, this was done by running the 
fireball an extra ~1 fm/c u sing the in-medium p spec- 



tral function. It turns out 19, Ml, however, that this 



description of p decays at thermal freezeout carries an 
extra factor of I/7 relative to a standard blast- wave spec- 
trum of hadrons at thermal freezeout, where 7 = qo/M 
is the usual Lorentz factor. Roughly speaking, the in- 
medium radiation given by Eq. (U) is proportional to the 
electromagnetic decay width times the fireball lifetime, 
(r(z/7)TFB, while decays after freezeout are proportional 
to the branching ratio of electromagnetic to total life- 
time, Tii/r^p, where the 7 factor cancels (as usual, we 
define Fx as the partial decay width in the rest system 
of the particle). The remainder of this section will give 
a more detailed discussion of this. 

To account for the correct time dilation effects in the 
calculation of dilepton decays of p-mesons after ther- 
mal freezeout, we use the standard Cooper-Frye descrip- 
tion Wdi, 



dN = q^da^' 



d'q 



(27r)3g' 



T 



(19) 



for the phase-space distribution of an on-shell particle 
(9^i'7^="^^) at thermal freezeout; T and = 7(1, /?i.) 
denote the local temperature and four velocity (flow) of 
the fluid cells of the medium, and da'' is the hypersur- 
face normal vector defined by an appropriate freezeout 
condition. In accordance with the homogeneous fireball 
model described in Sec. IIVI thermal freezeout at a con- 
stant time t = x'^ in the laboratory frame is assumed, i.e., 
(dcr^) = {d^x, 0, 0, 0). The in-medium p spectral function 
at freezeout is introduced via the substitution 



d^ 

90 



d^q 2S+{qf,,q^' -m^) 



j4 

d q — - 



(20) 



ip — -2/3 Im(_Dp)J^ (as before) includes the av- 



where Ap 

erage over the three polarizations via the Lorentz trace 
which runs over the physical {i.e., four-momentum trans- 
verse) components of the propagator. To properly treat 
the low-mass tails of the spectral function, one would 
have to resolve the individual resonance decays figuring 
into the p selfenergy. To simplify our task, we circum- 
vent this problem by employing the vacuum form of the 
p selfenergy augmented with a width corresponding to 
the full-width-half-maximum of the in-medium spectral 
function at thermal freezeout (F^" ~ 260 MeV). The dis- 
tribution of p mesons at thermal freezeout then reads 



d7Vf° 
d^xd'^p 



P 



T 



(21) 



In accordance with the averaging procedure over the 
freezeout duration, we evaluate the transverse four- 
velocity for the freezeout p at a time ATp°/2 ~ 0.5 fm/c 
after the thermal emission, Eq.([T]), has shut off. The 
dilepton rate follows by folding with the appropriate par- 
tial decay width. Within our vector-meson dominance 
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model this is given by the matrix element for the process 
p — > 7* + and after integration over t € (tfo, oo) 
this results in 



dN, 



(fo) 



T 



1 

pfo 

p 



go t 

Mrf° 



(22) 



where L is the dilepton-phase space factor ([3]). The 
second line of Eq. shows that the momentum de- 
pendence of the dilepton distribution from p decays af- 
ter thermal freeze-out deviates from the rate from a 
thermal source (H]) by a Lorentz factor 7 = qo/M ^ 
y^M^ + q^ /M . The physical origin of this difference is 
the time dilation of the total lifetime of a freeze-out p 
meson with three-momentum q which is absent in the 
formula for radiation from a thermal source, because its 
p-meson abundance at each instant of time is fixed by the 
temperature and pion-chemical potential of the medium 
(as required by detailed balance of p formation and de- 
cay), and thus the total number of thermal dileptons is 
determined by the lifetime of the fireball. Note, however, 
that the thermal rate is proportional to Tu with an asso- 
ciated time dilation factor I/7. The freezeout-p dilepton 
spectra are thus equivalent to standard blast-wave de- 
scriptions of stable hadrons [sst - 

Our default assumption for the radial profile of the flow 
field is a linear dependence on the radius according to 



(23) 



where R{t) ~ ro+a±P /2 is the radius of the fire-cylinder, 
= a±t its surface speed, and r < R{t) the radial 
coordinate of the fluid cell related to the volume ele- 
ment d'^a; ~ 2'Krdrdz (since we neither address azimuthal 
asymmetries nor peripheral collisions in this work, we ap- 
proximate the cross sectional area of the fire cylinder as 
a circle). 



B. Primordial p Mesons 

Another source of non-thermal dileptons is the de- 
cay of p mesons which originate from primordial hard- 
scattering processes and traverse the interaction zone 
without equilibrating. We evaluate this contribution 
within a schematic jet-quenching model (20l |. as follows. 

First, we construct the qx spectrum of primordial p 
mesons assuming a power law. 



1 dA^pri, 
qr dqr 



A 



{i+Bq^r 



(24) 



with parameters estimated from p-p scattering data [7l| , 
B = 0.525GeV~^ and a — 5.5. The total number of pri- 
mordial p's in A-A collision is determined based on the 



empirical freezeout systematics of light hadron produc- 
tion, i.e., we calculate the expected total number of p 
mesons in In-In at Tc ~ 175 MeV and match the norm 
of the primordial spectrum, Eq. ([M)) . to it. At the same 
time, we have to take care of the correct scaling proper- 
ties of the spectrum at high momentum: while the (total) 
yield at low qT is proportional to the number of partic- 
ipant nucleons, the high-f/T yield should scale with the 
number of primordial A^-A^ collisions. We implement this 
transition by a continuous linear switching between the 
two regimes over the range 1 GeV < qr < i GeV. 

Second, we implement a Cronin effect for A-A collisions 
by a "Gaussian smearing" of Eq. ((24l) , 



dA^"'°" 

pnm 

i2 

d qr 



d^kr dN, 



prim 



7rAfc2 



d^kn 



exp 



[qr - kxY 



Ak^ 



(25) 



with a conservative estimate of Ak^ =0.2 GeV^ as ex- 
tracted in Ref. [ll| based on direct photon spectra in p-A 
reactions. 

Third, we calculate a suppression factor representing 
the probability for primordial p mesons to escape the 
medium without rescattering, using Monte Carlo tech- 
niques. In line with our fireball model, we start from 
a spatially homogeneous distribution in the transverse 
plane at the QGP formation time. The escape probabil- 
ity for a p with momentum q is then calculated as 



P 



exp 



dt af^it) g{t) 



where the absorption cross section 



^abs 



Cph ~ 0.4 mb 
(Thad = 5 mb 



for 
for 



t < qo/mpTf 
t > qo/mpT[ 



(26) 



(27) 



depends on the p-meson formation time in its rest frame, 
which is assumed to be Tf=l fm, augmented by Lorentz 
time dilation in the fireball frame; (Tph and (Thad denote 
the pre-hadronic and hadronic absorption cross sections 
of the p, respectively, and g(t) the corresponding particle 
density of the fireball medium (partonic or hadronic; if 
the p meson has not formed by the time the system has 
hadronized, we use constituent quark scaling to infer the 
partonic density as = 3phad, and vice versa). 

A rough check of our construction may be obtained by 
comparing the p spectra to measured pion spectra in A-A 
collisions at the SPS. Since pion spectra contain both a 
hard (primordial) and a soft (thermal -t- flow) component, 
a comparison with our p spectra should contain both the 
primordial part constructed in this section and the ther- 
mal freezeout contribution described in the previous sec- 
tion, nil Al We do this comparison in terms of the nuclear 
modification factor, Raa, defined as the ratio of the spec- 
trum in A-A collisions over the collision-number scaled 
spectrum in p-p. The upper panel of Fig.[6]shows that the 
jet-quenching -|- freezeout-p model results in fair agree- 
ment with a recent compilation of Raa for pions in S-Au 
and Pb-Pb systems from various experiments at full SPS 
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FIG. 6: (Color online) Nuclear modification factor for trans- 
verse momentum spectra of p (upper panel) and n mesons 
(lower panel) from primordial (hard) production plus ther- 
mal freezeout, compared to a compilation of data for Raa of 
pions at full SPS energy [73 | for central collisions with the 
standard equation of state (EoS-A, full line) and an EoS with 
Tch = 160 MeV (EoS-B+C as discussed in Sec. |IV1 dashed 
line) . The primordial p spectra include effects of initial Cronin 
smearing and jet-quenching in the expanding fireball, while 
the thermal freezeout p's are taken from Sec. IIII Al (same for 
the pions, but without finite- width effects). 



energy [T^l- As a check of this procedure, as well as of 
the fireball model, we plot in the lower panel of Fig.[n]thc 
calculated pion Raa for central In-In at SPS (as a sum of 
the jet-quenching plus thcrmal-freezeout components) for 
central In-In at SPS; the agreement with data is reason- 
able (note that resonance decays are not included, which 
are mostly concentrated at low qx). 



depend on the strong coupling constant (it is purely elec- 
tromagnetic, (9(a°a^j„)) and can therefore be rather re- 
liably evaluated, provided one has a good knowledge of 
the parton-distribution functions (PDFs) within the nu- 
cleon. In a central collision {b = 0) of two equal nu- 
clei with mass number A, the invariant-mass spectrum 
of DrcU-Yan pairs per unit rapidity is given by 



dMdy 



fc=0 



47ri?2 



dAfdy 



(28) 



in terms of the standard DY cross section in an elemen- 
tary nucleon-nucleon (N-N) collision. 



^'' g^ [g(a;i)g(a;2) + g(a;i)q(a;2)] 



dMdy 

q—u,d.s 

(29) 

Here, q{xi^2) and q{xi^2) denote the (coUincar) quark 
and anti-quark distribution functions, respectively (ne- 
glecting nuclear effects). Their arguments are related to 
the center-of-mass (cms) rapidity, y, and the invariant 
mass of the lepton pair as xi^2 = xe^"^ with x = Mj 
where is the cms energy of the N-N collision. The 
root-mean-squared radius parameter in Eq. (|28p . i?o — 
1.05 fm, arises from a folding over a Gaussian thickness 
function; for simplicity, we will adopt Eq. (j28p also for 
non-central collisions with an accordingly reduced num- 
ber of participants, A. For the parton distribution func- 
tions we employ the leading-order parameterization from 
Ref. [zl] (GRV94LO), which incorporate isospin asym- 
metries in the sea-quark distributions which significantly 
reduce the need for additional isospin corrections for nu- 
clei with N ^ Z {e.g., less than 5% for Pb-Pb collisions), 
which will be neglected here (tJ- In Eq. higher- 
order corrections in as are encoded in an empirical K 
factor, which turns out to be iC ~ 1.5 to reproduce DY 
production in p-A collisions [zl] . 

In addition to generating a K factor, higher-order ef- 
fects manifest themselves in a nonzero pair-qx of the DY 
dileptons. To obtain a realistic qr spectrum, we follow 
the procedure adopted by the NA50 collaboration 0, [1| : 
based on a comprehensive analysis of p-A and A- A col- 
lisions at the SPS it has been found that the IMR and 
HMR dilepton spectra can be fairly well described by 
assuming a Gaussian shape for the DY spectrum. 



dTV, 



DY 



dA^DY CXp(-g2/2cr2 ) 



dMdydq:^ dMdy 



2< 



(30) 



C. Drell-Yan Annihilation and Correlated Charm 
Decays 

In the high-mass region (HMR, M > 3 GeV), dilep- 
ton spectra in nuclear collisions are expected to be dom- 
inated by the Drell-Yan (DY) process, i.e., primordial 
annihilation of quarks and antiquarks within the incom- 
ing nuclcons. To leading order, this process does not 



with 



'<3r 



0.8-1 GeV. In Ref. 



a value of a, 



9t 



0.8 GeV has resulted in reasonable agreement with the 
NA50 qr spectra for IMR dimuons in central Pb-Pb; for 
CTqy ~ 0.9 GeV, which increases the spectrum at qr = 
2 GeV by 50%, the agreement with NA50 would stiU be 



acceptable. Here, we adopt a„ 



0.8 GeV providing a 



conservative estimate of the DY contribution. 

The extrapolation of the DY spectrum to the regions 
where both mass and momentum are small (say, below 
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M,qT — 1.5 GeV) is problematic, but fortunately its 
contribution in this regime is small compared to, e.g., 
thermal emission in A-A collisions. There is, however, 
an additional constraint provided by the photon point 
(M 0) which allows to extrapolate the DY spec- 
tra to low mass, at least for reasonably large momenta, 
gT > 1 GeV (which is the regime were the DY yield is no- 
ticeable). Formally, the photon-production rate follows 
from the dilepton one by taking the limit M — > 0. More 
specifically, this is encoded in the relation 



9o 



d^ 



^— Imn(^„,,^(M = 0,q)f^iqo;T) (31) 



with the same e.m. current-current correlation function, 
n^j'^j of Eq. (HI), as in the dilepton rate, Eq. (P). For 
q ^ M, the M-dependence of Ilom is weak, and the dif- 
ference between qo ^^j^"^ and go amounts to a factor 

of Thus we can evaluate the DY qx spectrum at a 

mass Mcut and extrapolate it down in mass to the photon 
point using the factor Mcut/M. For Mc„t=0.8-1 GeV, in 
connection with (jq^ ~ 0.8 GeV, reasonable agreement 
with the primordial photon spectrum of Ref. [ll| is found. 

In addition, there could be "pre-equilibrium" contri- 
butions from secondary Drell-Yan processes [zl] {e.g., 
■nN jjL^X involving primordially produced pions), 
which turn out to be rather sensitive to the pion for- 
mation time. E.g., for t^^^^ = 1 fm/c, the enhance- 
ment over primordial Drell-Yan annihilation in central 
S-U was found to be -10% at M = 2 GeV. In a ther- 
mal emission description, which we employ here after a 
rather early thermalization time of tq = 1 fm/c, it is dif- 
ficult to separate pre-equilibrium radiation from thermal 
emission (e.g., secondary Drell-Yan in 7r-iV interactions 
might overlap with thermal q-q annihilation). A rough 
estimate of pre-equilibrium effects may be obtained by 
varying the thermalization time; decreasing, e.g., tq to 
0.8 fm/c (which is close to the overlap time of the 2 col- 
liding nuclei at SPS energies) increases the QGP contri- 
bution at = 2 GeV by ~50%, which is less than 10% 
of the Drell-Yan contibution at all masses considered. 

For the dilepton contribution from correlated decays 
of D and D mesons we use the experimental result 
from p-p collisions extrapolated to In-In, as provided 
by the NA60 collaboration [l^. As a note of caution, 
we remark that recent measurements at the Relativistic 
Heavy-Ion Collider (RHIC) report substantial modifica- 
tions of heavy-quark spectra in Au-Au collisions, rela- 
tive to p-p (as inferred from a suppression and elliptic 
flow of "non-photonic" single-electron spectra associated 
with semileptonic decays of open-charm (and -bottom) 
hadrons) [t^, [t^. Such (possibly nonpcrturbative fl% ) 
medium modifications of the charm momentum spectra 
presumably translate into a softening of the invariant- 
mass spectra of pairs as well. At the SPS, the 
shorter QGP lifetime is likely to lead to smaller effects, 
but an explicit measurement of the (delayed) charm de- 
cays has been presented recently [ts} . 



IV. THERMAL FIREBALL EVOLUTION 

Our description of the space-time evolution of central 
and semicentral A-A collisions is approximated by an ex- 
panding thermal fireball characterized by a time depen- 
dent cylindrical volume as [23 |. 



1 



1 . 

^0 + Vz.nt + -jayV 



(32) 

where we neglect effects due to a finite ellipticity. The 
initial transverse radius r_L^o is determined by the central- 
ity of the collision {e.g., rj^.o = 5.15(4.6) fm for central 
(semicentral) In-In collisions). The initial longitudinal 
size, zq, is equivalent to the formation time, tq, of the 
thermal medium, which we fix at the standard value of 
To = 1 fm/c (translating into zq ~ To Ay = 1.8 fm, where 
Ay = 1.8 represents the rapidity width of a thermal fire- 
ball). For the longitudinal expansion we employ a mod- 
erate acceleration, — 0.045c^/fm, together with an 
initial velocity of w^.o = 0.6c (reminiscent to Ref. [2^). 
but the dilepton invariant-mass spectra are essentially 
unaffected if we use [vz^^ii, fflz) = (c, 0) as in Ref. [11]. The 
most important parameter is the transverse acceleration. 
More recent applications of the fireball model, both in the 
dilepton [lH| and charm diffusion IttI c ontext, have used 
larger values than in previous work [24j. [38l | , in the range 
aj^ = 0.08 — O.lc^/fm. Here, we employ a\_ = 0.085c'^/fm 
as in Ref. 

The time evolution of the temperature is determined 
assuming entropy conservation. At a given collision en- 
ergy, the hadro-chemistry of the fireball is characterized 
by thermal- model fits to the observed hadron ratios. In 
our default scenario we assume the chemical freezeout 
temperature to coincide with the critical temperature for 
QGP formation, at Tch = = 175 MeV. For central 
and semicentral In(158A GeV)-|-In collisions we fix the 
entropy per (net) baryon at s/g^*' = 27, which, using 
a hadronic resonance gas (HG) equation of state (EoS), 
translates into an associated baryon chemical potential 
of iJ.'j^ = 232 MeV, well within the uncertainties of re- 
cent thermal model fits at SPS [2l|, [2^ (all other chem- 
ical potentials being zero). The subsequent hadronic 
trajectory in the phase diagram is then constructed at 
fixed s/g^g^ = 27 (isentropic expansion) with the addi- 
tional constraints of pion, kaon, 77 and antibaryon num- 
ber conservation, requiring the build-up of corresponding 
chemical potentials. For a given number of nucleon par- 
ticipants [i.e., collision centrality), the fireball entropy 
amounts to = 2630(1890) for central (semicentral) In- 
In collisions, translating into a charged particle multi- 
plicity of dNch/dy ~ 195(140) ~ A'part (for a chemical 
freezeout temperature of Tch = 160 MeV, as considered 
below, the multiphcities increase by less than 5%). Us- 
ing the relation between the total entropy and volume, 
S = s{t) ■ V{t), the entropy density, s{t), can be used to 
determine the time evolution of temperature and baryon 
density (along the hadronic trajectory) using the HG 
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EoS, and similarly for the QGP phase using a quasipar- 
ticle EoS. The volume partition in the HG-QGP mixed 
phase is calculated from the standard mixed-phase con- 
struction, where the fraction of matter in the hadronic 
phase is given by 



.QGP 



sit) 



oQGP 



„HG 



(33) 



where sHG.qgp denote the critical entropy densities at 
Tc- With a formation time of tq = 1 fm/c, the evolution 
for central (semicentral) In-In starts in the QGP at To = 
197(190) MeV, passes through a mixed phase at T'^^ = 
Tc = 175 McV, and terminates at thermal frcezcout at 
around Tfo ~ 120-135 MeV. 

The main uncertainties associated with the fireball evo- 
lution are the transverse acceleration as well as the overall 
lifetime (which is somewhat correlated to the longitudinal 
expansion). It turns out that the latter is rather sensi- 
tive to the absolute magnitude of the experimentally ob- 
served dilepton excess radiation, resulting in Tfo ~ 6 fm/c 
for central In-In collisions. The remaining uncertainty 
consists of the interplay between longitudinal and trans- 
verse expansion. As mentioned above, the variations 
considered in our previous works [H, [ISl, a_L = 0.08- 
0.085 cVfm and {v°,a^) = (c, 0)-(0.6c, 0.045 c^/fm), have 
negligible impact on the invariant mass spectra, while 
the longitudinal-acceleration scenario, which we will fo- 
cus on here, allows for a slightly larger transverse ex- 
pansion which appears to be favored by the rather hard 
qx spectra. In principle, more accurate information on 
the final state of the expansion can be obtained once 
hadronic spectra for In-In are available (our pion spec- 
tra are actually in reasonable agreement with S-Au data, 
recall lower panel of Fig. [6]). For central Pb-Pb colli- 
sions at SPS energies, the typical results for transverse 
surface velocity and thermal freezeout temperature are 
(/31/c,Tfo[MeV]) = (0.65±0.1,120±10) [H, H [Ml, HI , 
which is quite comparable to the values used here for In- 
In. We emphasize that all contributions to the dilepton 
spectrum (QGP, p, uj, (p, and four-pion) are tied to the 
same evolution, thus fixing their relative weights. 

Finally, to illustrate the uncertainties associated with 
the underlying equation of state we will investigate three 
combinations of critical and chemical freezeout tempera- 
tures roughly covering the current theoretical and exper- 
imental ranges: 

(A) Our default scenario, employed in most of our cal- 
culations thus far, consisting of identical Tc and 
Tch at an "intermediate" value of 175 McV; ther- 
mal freezeout is fixed at (/Lt^", Tfo) = (79, 120) MeV 
(semicentral In-In). 

(B) A scenario with a relatively small and identi- 
cal critical and chemical freezeout temperature at 
(Mb > ^ch) = (240, 160) MeV, compatible with re- 
cent thermal model fits in Refs. [Ill, [22|. For 
dilepton spectra the most important consequences 



of this scenario are a significantly extended QGP 
phase which will increase its thermal emission con- 
tribution and reduce the hadronic one, in partic- 
ular at intermediate masses. In addition, due to 
smaller pion chemical potentials in the subsequent 
hadronic phase, the freezeout temperature (at fixed 
fireball lifetime) will be larger than with EoS-A, at 
(Mi°,Tfo) = (37,136) MeV. 

(C) A scenario with a large critical temperature Tc = 
190 MeV (as sug gested by recent lattice QCD 



computations |23j . which maximizes (minimizes) 



the space-time volume occupied by the hadronic 
(QGP) phase. Since chemical freezeout at such a 
temperature is questionable, we allow for a chem- 
ically equilibrated hadronic phase until chemical- 
freeze out sets in under the same conditions as in 
EOS-B, at (/if, Tch) = (240, 160) MeV. 



V. COMPARISON TO DILEPTON SPECTRA 
AT SPS 

We now turn to a systematic analysis of experimental 
dilepton spectra as measured at the SPS by the NA60 
and CERES/NA45 collaborations in In-In and Pb-Au 
collisions, respectively. Based on the various ingredients 
developed in the previous sections, we first address the 
NA60 invariant-mass and transverse-momentum spectra, 
followed by a consistency check with earlier and updated 
CERES data. 



A. Invariant Mass Spectra 

Thermal pT invariant-mass spectra for A- A col- 
lisions are computed by integrating the emission rate, 
Eq. ([1]), over the fireball evolution (as well as three- 
momentum). 



dM Ay J J qa d^x d^q 



(34) 

where Acc denotes the detector acceptance which has 
been carefully tuned to NA60 simulations [s^. The 
fugacity factor, Zp=e"''^/^, arises due to chemical off- 
equilibrium in the hadronic phase for T < Tch; it de- 
pends on the thermal emission source under considera- 
tion, cf. Eq. (I4|): for the p, to and four-pion contribu- 
tions one has with n = 2,3,4, respectively, while 
for the (f) one has -7^, where 7^ ~ 0.75 accounts for 
strangeness undersaturation for medium-size nuclear col- 
lision systems at the SPS [sj. In addition, appropriate 
fugacities figure into the various in-medium sclfenergy 
contributions. 

Initial comparisons of NA60 data to theoretical 
predictions |14| have focused on the contribution from 
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FIG. 7: (Color online) NA60 excess dimuon spectra [T^l in 
semicentral (upper panel) and central (lower panel) In-In col- 
lisions at SPS compared to theoretical calculations using an 
in-medium e.m. spectral function. The individual contribu- 
tions arise from in-medium p-mesons [SJl (dash-dotted red 
line), 4n annihilation with chiral V-A mixing (dashed blue 
line) , QGP plus correlated open charm decays (dotted orange 
line) and Drell- Yan annihilation (solid turquoise line) ; the up- 
per dashed brown line is the sum of the above, while the solid 
purple line additionally includes in-medium ui and (p decays as 
well as freezeout plus primordial ps (solid black line). In the 
semi-central data, the uncertainty due to the -q cocktail sub- 
traction is indicated by the open and filled data points (the 
former are based on an estimated rj yield at high qr while the 
latter represent an upper limit on the 77 by subtracting the 
dimuon spectrum to zero at threshold). 



the p meson which dominates in the LMR. The shape of 
the in-medium p-spectral function describes the experi- 
mental spectra well, but the absolute yields have been 
overestimated by ~ 30%. This discrepancy has been re- 
solved [iBl by increasing the transverse fireball expansion 
{a±), reducing the fireball lifetime to about 6-7 fm/c, 
cf. Sec. HV] above. In addition, the larger transverse ex- 
pansion leads to harder emission spectra in qt, which will 
be helpful in understanding the qx spectra, as discussed 
in the following sections. 

Fig. [7] summarizes our results for the mass spectra in 
semicentral and central In(158 AGeV)-In collisions, com- 



FIG. 8: (Color online) Upper panel: NA60 data [T3| com- 
pared to thermal dimuon spectra using (i) in-medium p-, 
uj- and (^-mesons without baryon effects (-|-QGP-|-charm-|-in- 
medium four-pion; solid purple line), and (ii) free p 
(-|-QGP-|-charm-|-free four-pion; dashed green line). Lower 
panel: NA60 data [l^ in the IMR compared to thermal 
dimuon spectra with different implementations of the four- 
pion contribution, using either its vacuum form (lower dashed 
line) or including chiral mixing (upper dashed line), and cor- 
responding total spectra (lower and upper solid line, respec- 
tively) . 



puted for the EoS-A scenario (Tc T^h = 175 MeV). The 
modifications relative to our previous work are: (i) 
the freezeout p has been separated from the in-medium p 
contribution, (ii) primordial p and (iii) Drell- Yan contri- 
butions have been added. As a result of the separation 
(i), the thermal emission lifetime is now slightly smaller, 
6.5(6.2) fm/c for central (semicentral) collisions. The in- 
clusive mass spectra (and pertinent conclusions) are es- 
sentially identical to those in Ref. [l^, i.e., the predicted 
in-medium effects on p spectral function lead to good 
agreement with the data in the LMR, while the same 
fireball evolution also reproduces the observed excess in 
the IMR well. The largest source here is four-pion anni- 
hilation, together with smaller contributions from open- 
charm decays and QGP emission. In-medium (f) decays 
are noticeable but not very significant relative to the un- 
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FIG. 9: (Color online) Dimuon invariant mass spectra in two 
different bins of transverse pair momentum, qr [13 ■ Upper 
panel: qr < 0.5 GeV; lower panel: > 1.0 GeV. 



FIG. 10: (Color online) The same as Fig. |9] but for central 
collisions. 



certainty in the data. 

The sensitivity to the medium effects in the e.m. cor- 
relator is further illustrated in Fig. [51 The upper panel 
demonstrates that a free p spectral function is ruled out, 
but also one which only includes modifications due to 
a meson gas clearly does not reproduce the data, due 
to both a too narrow peak and a lack of enhancement 
below the free p mass, especially when approaching the 
dimuon threshold. In the lower panel one sees that the 
effect of chiral mixing on the 47r contribution amounts 
to up to a factor of ~2 enhancement in the ai resonance 
region, but the effect on the total is rather moderate. 
Thus, no strong case on the chiral mixing can be made 
at present. The prevalence of the four-pion contribution 
in the IMR is reminiscent to the hydrodynamic calcu- 
lations in Ref. [lB| where hadronic rates calculated in 
the chiral reduction formalism have been employed. On 
the contrary, in the fireball calculations of Ref. [13| QGP 
emission dominates in the IMR. We elucidate on this dis- 
crepancy in Sec. IV CI below. 

In Figs, ini and [TO] we compare our theoretical calcu- 
lations to NA60 Af-spectra binned into regions of low 
{qt < 0.5 GeV) and high {qt > 1.0 GeV) transverse pair 



momentum, for both semicentral and central collisions. 
Also in this representation the agreement is fair. There 
is possibly an indication of a slight over- (under-) esti- 
mate in the high-qy bin for (semi-) central collisions. To 
further scrutinize this issue we now turn to qt spectra, 
binned in invariant mass. 



B. Transverse Momentum Spectra 

In analogy to Eq. ([M)) for AI spectra, qx spectra are 
computed by integrating the dilepton rate, Eq. ([1]), over 
the space-time evolution of the fireball and a suitable in- 
terval in invariant mass, [Mmin, il/max]- However, since 
the available experimental spectra have been corrected 
for the detector acceptance, the Monte-Carlo simula- 
tion of the acceptance function (with radial and angular 
dependencies due to directed (radial flow) and random 
(thermal) motion of the virtual photon) is now replaced 
by explicit integrations over the radial coordinate, r, of 
the fireball, the azimuthal angle, c/)q, of the virtual photon 
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momentum and its rapidity, y, 



In the non-relativistic limit, qr <C M, we use the asymp- 
totic form Iq{x) =x-tO 1 + ~ exp(x^/4) to obtain 



dNu 



tto RW M„ 



2iT y> 



J dt J dr j dM j d(j)q J dy2TTrz{t) 



M„ 



y< 



xM z 



dNu 



P j4 i4 

a xa q 



(35) 



As before, R{t) denotes the radius of the expanding 
fire-cyhnder (corresponding to the first parenthesis in 
Eq. ([5^ ). z{t) its longitudinal length (second parenthe- 
sis in Eq. ([5^ ). Note that the rate, Eq. is calcu- 
lated in the thermal rest frame while the integrations in 
Eq. ((35)) are in the laboratory frame. The relation be- 
tween the four-momentum, g^, in the laboratory frame 
and the one in the local rest frame of the fluid cells, g^, 
is determined by a boost with the radial flow velocity, 
Eq. ((23|) . Due to rotational invariance the thermal rate, 
Eq. Il]), depends only on the magnitude of the three- 
momentum in the thermal rest frame and therefore one 
of the integrations over the spatial and momentum az- 
imuthal angles in Eq. (j35[) becomes trivial yielding a fac- 
tor of 27r. To gain qualitative insights into the behavior of 
dilepton qT spectra from thermal sources, let us assume 
that the three-momentum dependence of the hadronic 
e.m. current correlator is weak (it would be absent in 
the absence of medium effects), so that for not too large 
invariant-mass intervals, the integrand in Eq. (|35p can 
be considered constant except for the Bose distribution. 
For M ^ T, it is possible to approximate the qualitative 
behavior of the qx spectra in analytic form [s^, [s^ . For 
a linear flow profile, Eq. (|23p . and in Boltzmann approx- 
imation one has 
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=C I dr t 
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(36) 



with C weakly qx dependent. Substituting = 
TOj_ cosh y {rax = \J A/^ + g|,) the integrals over y and 
are given by modified Besscl functions. For the y-integral 
we can use their asymptotic form for vtit jT > AI/T ^ 1: 
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X exp 



(38) 



We have numerically verified that this expression can 
be further simplified by the following approximate treat- 
ment of the r integral: In the pre-factor of the exponen- 
tial substitute /3± by its average over r, 



2 

R^ 



drrPi^^^Pi; (39) 



use the r average of the argument in the exponential, 
which, together with the non-relativistic approximation 
qrp =g^^]^,j = yj2M{mT — A/), results in 



qjqrdMdt ,.|M^r 



X exp 



272 



(40) 



where 



T'nff — 



1 / TTlT 

X \ I exp 

rriT \ Toff 



l-(M-T)(/31)2/(4r) 



(41) 



and {(3\)^ = (/31)V2. The form for /3i < 1 is the 
known "pocket formula" for the parametric dependence 
of the effective slopes of hadronic qx spectra on the fire- 
ball temperature, the particle mass and the surface-flow 
velocity [8^|. We recall that in the ultrarelativistic limit, 
gT 3> Af, the dependence of the effective temperature 
on the radial velocity is given by the Doppler blue-shift 
expression for a thermalized gas of massless particles. 



Teff(/3±) 




(42) 



which follows from the large-gr limit of /q in Eq. ([57]) [86j 
which also provides an additional factor \/ ^/qr. For our 
r-dependent flow profile the effective temperature is given 
by the blue-shift value for an average value (/3j_)^ = 
where ^ G (0, 1) is not easily determined by simple ap- 
proximations of the radial integral. Numerical stud- 
ies show that in the case of a linear radial-flow profile, 
Eq. (|23p . and typical parameters in the region of the p 
peak (Af = 0.75 GeV, (31 ~ 0.5), i ~ 0.8-0.85 leads to a 
good estimate for T^s at high qr- 
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It is important to note that for vector-meson decays 
after thermal freezeout, of. Sec. IIII Ap . the qt spectra 
are harder by an additional factor ttit/M, i.e., applying 
the same approximations to Eq. (j22p as used to derive 
Eq. ((40l) . one obtains 



dN, 
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X exp 



M2(/31)2 



2r2 



(43) 



M 
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nix 
Tcff 



.(fo) 



with Toff again given by Eq. ([iT|) . 

The additional factor tot /M in ((43|) compared to ([40]) 
originates from the additional factor jM in the Cooper- 
Frye (CP) formula ([22)) compared to the McLerran- 
Toimela (MT) emission formula ([1]), as discussed in detail 
in Sec. [Tira 

The above discussion of effective slopes in thermal 
emission spectra obviously does not apply to non-thermal 
dilepton sources, i.e., decays of "primordial p mesons" 
fSec lmBj) and Drell-Yan annihilation fSec HlTC)) . Both 
are characterized by a power-law behavior at high qy, 
where their contribution becomes potentially important. 
Toward lower ^t, the primordial p contribution is much 
suppressed due to p absorption ("jet quenching"), while 
the Drell-Yan process is no longer well-defined. How- 
ever, the Drell-Yan dileptons carry the hardest slope of 
all sources considered, which renders even a naive ex- 
trapolation of their spectra to low qr very small (no more 
than a few percent for < 1 GeV in all mass bins below 
M = 1.4 GeV). M QT ~2 GeV, in turn, the Drell-Yan 
contribution is quite appreciable. 

Figs. [TT] and [T^] summarize our calculations for semi- 
central and central In-In collisions, respectively, in three 
different mass bins in comparison to the NA60 data [3 ■ 
Open-charm decays have been removed from the exper- 
imental spectra and are consequently not included in 
the theory curves either. The parameters of our fire- 
ball evolution are as described in Sec. IIVI within our 
default EoS-A scenario (T^ = Tch = 175 MeV). Since 
the overall normalization of the experimental spectra is 
not known at present, we have fixed it in each mass 
bin using the invariant mass spectra in the low-g^ bin 
(0 < (7t < 0.5 GeV). This easily translates into a 10- 
20% uncertainty in the absolute normalization. For the 
central collisions the agreement between theory and data 
is quite satisfactory in all mass bins. In the high-mass bin 
1 GeV<M<1.4 GeV, lower panel in Fig. [H]), this pro- 
cedure appears to overestimate somewhat the gr spectra 
aX qx >1 GeV, which, however, seems not to be refiected 
in the corresponding qx bin in the mass spectra (lower 
panel in Fig. [TU)) . This could very well be due to the 
uncertainty in the underlying normalization procedure. 
Note again that the compared to Ref. [l^ additionally 
implemented hard components (primordial p's and Drell- 
Yan) arc insignificant at < 1 GeV, i.e., for the total 
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FIG. 11: (Color online) Dimuon transverse pair-momentum 
spectra in semicentral In(158 AGeV)-In collisions for three 
bins of invariant mass [lg|. Upper panel: 0.4 GeV < M < 
0.6 GeV; middle panel: 0.6 GeV < M < 0.9 GeV; lower panel: 
1.0 GeV < M < 1.4 GeV. 



yield and the understanding of the spectral shape of the 
inclusive A'/-spectra. On the other hand, at higher qx, 
these contributions are essential (even dominant) for a 
proper description of the spectra. In fact, in semicentral 
collisions, the experimental transverse momentum spec- 
tra in the two mass bins below (0.4 GeV< qx < 0.6 GeV) 
and around (0.6 GeV< qx < 0.9 GeV) the free p-mcson 
mass turn out to be underestimated (and too soft) for 
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(Color online) Same as Fig. 1111 but additionally 
(factor 2 augmented) contributions from uj f-channel 
in Tvp ■Rfip. 



> 1 GcV. Especially in the low-mass bin, the discrep- 
ancy again appears to be larger than one could anticipate 
from the mass spectrum in the 0.4-0.6 GeV regime in the 
(7t > 1 GeV bin (lower panel in Fig. [9l). The fact re- 
mains, however, that the theoretical qx spectra in the 
low and intermediate mass bin are somewhat too soft. 
Even though the final transverse flow velocity of the fire- 
ball model is about 7% smaller for semicentral relative to 
central In-In collisions, this difference would not be able 
to account for the discrepancy (it amounts to a change 
in slope by about 10 MeV for the freeze-out p, and even 



less for thermal radiation). The nuclear suppression fac- 
tor, RAA{qT)i in Fig. [B] (solid hue) suggests that there 
might be room for a 30% increase of the primordial plus 
freezeout p contribution (possibly more if one could com- 
pare to the Raa of P mesons; due to their larger mass a 
more pronounced flow effect could enhance their Raa)- 
Whether this would suffice in the low-mass bin, is ques- 
tionable. 

Another possibility consists of meson t-channel ex- 
change in elastic scattering of p mesons off pions in 
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the heat bath (with subsequent conversion into a Icp- 
ton pair). These processes have been calculated previ- 
ously for real photon production [ll[ (which arc given 
by the same diagrams except that the final-state pho- 
ton is on-shell), where exchange has been found to be 
the most important process at high qr', however, con- 
tributions from other exchange processes (e.;;., tt and 
fli exchange in Trp ^ 7r7, or reactions involving strange 
mesons [ll[) are also appreciable. To roughly account for 
the latter, we augment the rates computed in Sec. Ill Dl 
by a factor of two (in addition, we recall that Trp initial 
states in the hadronic fireball emission carry a pion fu- 
gacity factor to the third power, zW The convolution 
of these rates over the fireball evolution leads to a con- 
tribution to dilcpton-gr spectra as shown by the lower 
dash-dotted line in each of the panels in Fig. [T31 it in- 
deed provides the hardest spectrum among all thermal 
sources (QGP, in-medium vector mesons and four-pion 
annihilation) , in line with the rates displayed in Fig. [S] 
Consequently, the relative importance of the ^-channel 
processes grows with qr but remains rather moderate 
even at — 2 GeV, up to 15% and 5% of the total the- 
oretical yield in the M = 0.4-0.6 GeV and 0.6-0.9 GcV 
mass bins, respectively (negligible for M > 1 GeV). Note 
that the slopes of the ^-channel emission spectra resem- 
ble the data quite well, but our present estimate of their 
strength is insufficient to resolve the discrepancies at high 
qT- However, their impact on the slope of the total spec- 
tra is not insignificant, as we will sec in Sec. IV E)l below. 



C. Hadro-Chemical Freezeout and Critical 
Temperature 

In all our calculations of thermal dilcpton spectra thus 
far, the medium evolution was based on the notion that 
the critical (Tc) and chemical freezeout (Teh) tempera- 
ture coincide, at = Teh = 175 MeV. More recent the- 
oretical (lattice QCD) and phenomenological (thermal 
model fits to hadron ratios) studies, however, allow for 
the possibilities that Tc could be significantly larger (190- 
200 MeV) and Teh significantly smaller (150-160 MeV). 
While a smaller — 160 MeV is still viable, a larger Tch 
appears unlikely since in a high-density hadronic phase 
number-changing reactions affecting the chemistry are 
to be expected (in addition, thermal model fits start 
to become unstable at temperatures above ~180 MeV 
due to uncertainties in the high-lying, high degeneracy 
hadron resonance spectrum). Therefore, in this sec- 
tion we study the sensitivity of the invariant-mass and 
qT spectra to the hadro-chemistry of the fireball, keep- 
ing its geometry and flow parameterization as well as 
the total lifetime the same as before (including all nor- 
malizations, where apphcable). As quoted in Sec IIVI 
in addition to our standard equation of state (EoS-A: 
Tc = Tch = 175 MeV), we investigate two alternative 
scenarios: EoS-B with = Tch = 160 MeV and EoS-C 
with Tc = 190 MeV, Tch = 160 MeV. In general, a lower 



chemical freezeout temperature (where all meson chemi- 
cal potentials are by definition zero) entails a larger vol- 
ume at chemical freezeout, and thus also smaller baryon 
chemical potentials (ie., smaller baryon densities), as 
well as smaller meson chemical potentials in the subse- 
quent evolution. This furthermore implies that, to obtain 
roughly the same overall dilcpton yield, thermal freeze- 
out occurs at a larger temperature (and smaller (Itt^k), 
which, in principle, leads to (somewhat) harder qx spec- 
tra, both at thermal freezeout and for thermal emission 
throughout the hadronic evolution. In the remainder of 
this Section, we focus on semicentral collisions. 

We start our investigation within the "EoS-B" sce- 
nario, i.e., Tc = Tch = 160 MeV. The rcsuhing M- and 
(7T-spectra are compiled in Figs. [Ml and [T5l respectively. 
The smaller Tc implies a longer duration of both the QGP 
and mixed phase (the latent heat for EoS-B is larger than 
for EoS-A since entropy density for the QGP EoS drops 
slower with T than for the HG EoS), while the duration 
of the hadronic phase is accordingly reduced. As noted 
above, the thermal-freezcout temperature has increased 
to Tfo = 136 McV (compared to Tfo = 120 MeV for EoS- 
A). As an immediate consequence, QGP radiation is sig- 
nificantly enhanced, while the hadronic yield is slightly 
reduced (both for in-medium p in the LMR and four-pion 
in the IMR). The overall quality in the description of the 
invariant-mass spectra is quite comparable to the EoS-A 
scenario. The somewhat harder spectra implied by the 
larger hadronic temperatures for EoS-B (most notably for 
the freezeout p) lead to a slight increase of the dilcpton 
yield in the p mass region in the qx > ^ GeV momentum 
bin (compare lower panels in Figs. and [H)) . However, 
in the low- (0.4 GeV< M < 0.6 GeV) and intermediate- 
mass (0.6 GeV< M < 0.9 GeV) bins of the qr spectra, 
the improvement at high qx is rather marginal, while in 
the higher mass (1.0 GeV< M < 1.4 GeV) there is es- 
sentially no change compared to the EoS-A calculation 
(recall Fig. [TT|) . despite the fact that the QGP contribu- 
tion is now larger than the four-pion one. 

Turning to the EoS-C scenario (summarized in Figs. 1161 
and[T71), the large value of = 190 MeV leads to a QGP 
and mixed phase which is shorter than for EoS-A, but 
there is now a "high-density" (chemically equilibrated) 
hadronic phase down to Tch = 160 MeV, followed by 
a subsequent chemical off-equilibrium evolution which is 
identical to EoS-B. Obviously, the QGP emission yield 
is the smallest of all three scenarios, while the hadronic 
yield (both in-medium p's and four-pion annihilation) is 
the largest, possibly implying a slight overestimate in the 
p-mass region of the inclusive mass spectra (13% larger 
than for EoS-A, which could be readjusted by a slightly 
reduced lifetime). The high-density hadronic phase helps 
in the qx spectra in the low- and intermediate-mass bin, 
for which, at high qx, the discrepancy is the smallest 
for the three scenarios. In the highest mass bin, even 
the rather large ratio of four-pion annihilation over QGP 
emission does not upset the agreement with the qx spec- 
tra, indicating that a large part of the four-pion emission 
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FIG. 14: (Color online) The same as Fig. [9l but for a fireball 
with a critical temperature Tc — 160 MeV and a chemical- 
freezeout temperature T^h = 160 MeV (EoS-B). 



FIG. 15: (Color online) The same as Fig.[TT] but for a fireball 
with a critical temperature Tc = 160 MeV and a chemical- 
freezeout temperature Teh = 160 MeV (EoS-B). 



emanates from the high-density hadronic phase. 

Since the question of four-pion vs. QGP emission in 
the IMR has drawn considerable attention in the recent 
literature [s^l, we take a closer look at their interplay 
in our three scenarios in Fig. [181 The plot illustrates 
again that with EoS-B (T^ = T^h = 160 MeV) the QGP 
contribution exceeds the in-medium four-pion annihila- 
tion over the entire IMR considered (note that the tails 
of the p spectral function are not included). The op- 
posite trend applies to our default EoS-A scenario, and 



even more so for EoS-C. Especially the large difference 
between EoS-B and EoS-C clearly demonstrates that the 
major portion of the four-pion emission emanates from 
the high-density hadronic phase (EoS-B and EoS-C have 
an identical hadronic evolution below Teh = 160 MeV); a 
similar argument applies to EoS-A, i.e., four-pion emis- 
sion does not primarily arise in the late stages (with large 
flow) despite an enhancement through large pion-fugacity 
factors. The basic reason for this can be understood 
by inspecting the temperature-mass differential emission 
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FIG. 16: (Color online) The same as Fig. [9l but for a fireball 
with a critical temperature Tc — 190 MeV and a chemical- 
freezeout temperature T^h = 160 MeV (EoS-C). 
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FIG. 17: (Color online) The same as Fig.[TT] but for a fireball 
with a critical temperature Tc = 190 MeV and a chemical- 
freezeout temperature Teh = 160 MeV (EoS-C). 



yield which is roughly given by [l^ 

(X IraU^M, T) e"^^/^ T'^-^ , (44) 

characterized by the standard thermal exponential fac- 
tor and a power in temperature resulting from the three- 
momentum integral over the Boltzmann factor and, most 
importantly, the volume expansion. Assuming a weak 
temperature dependence of the e.m. spectral function 
(medium effects bias the emission toward higher T), a 



differentiation over T identifies the temperature of max- 
imum emission as T,nax — M/5.5. In the IMR of interest 
here, this means a T^ax right around (note that, at 
Tc, the hadronic phase has an extra "advantage" of a 
roughly twice larger volume compared to the QGP due 
to the latent heat), consistent with our fireball emission 
pattern. 

We conclude this section by noting that "reasonable" 
variations in hadro-chemical frcezeout and critical tem- 
perature have a very moderate impact on the invariant- 



23 



semicentral In-In 
all p-j- 



QGP (EoS-A) 
QGP (EoS-B) 
QGP (EoS-C) 
47t (EoS-A) 
4it (EoS-B) 
4it (EoS-C) 



0.9 



1.1 1.2 

M (GeV) 



1.3 



1.4 



10 



10"^ 



model 

slope fit (0.4 GeV<qT,<1.8 GeV) 
slope fit (1.0 GeV<q.p<1.8 GeV) 




0.5 



1.5 



qx (GeV) 



FIG. 18: (Color online) Comparison of the QGP and four-pion 
contributions to the dilepton yield in the intermediate mass 
region for the three equations of state, as discussed in the text. 
The 2-pion part, which is significant up to M ~ 1.1 GeV, is 
not shown. 



FIG. 19: (Color online) Dependence of the effective-slope fits, 
Eq. (|45|l . to our theoretical qr spectra without the Drell-Yan 
contribution (dots) on the fit range in qr (as indicated in the 
legends) for the mass bin 0.3 GeV <M< 0.4 GeV. 



mass spectra as seen by NA60. On the one hand, this 
implies little sensitivity to the concrete values of and 
Teh- On the other hand, it means that our results are very 
robust with respect to uncertainties in these quantities. 
The main reason for this robustness is the "quark-hadron 
duality" of the underlying (medium-modified) emission 
rates from the hadronic and QGP phase in the relevant 
temperature regime, T = 160-190 MeV, in both LMR (p 
melting) and IMR (including chiral mixing). As such, it 
provides additional support to the medium modifications 
in the employed rates. In the IMR, this "duality" does 
not allow for a (maybe even academic) distinction be- 
tween a high-density hadronic or partonic source. Parti- 
tions with either component dominant are viable in both 
mass and transverse momentum spectra. Concerning the 
high-QT region of the lower two mass bins, our calcula- 
tions indicate a slight preference for EOS-B and EOS-C, 
due to their smaller Teh, implying larger temperatures in 
the hadronic evolution. 



D. Slope analysis 

To complete our analysis of the qx spectra we perform 
in this section a quantitative analysis of effective slope 
parameters following the procedure adopted by the NA60 
collaboration for the experimental data [l^, [83, [111 . The 
qt spectra have been divided into several invariant-mass 
bins, in each of which the data have been fitted to the 
function 

1 diV 1 dA^ / mT\ .... 

— — = = Cexp -— , (45) 

qt aqr mr drriT \ IcsJ 

where the fit range is typically taken as 0.4 GeV < qt < 
1.8 GeV, but the extracted slopes are rather insensitive 
to variations as long as the lower limit is above 0.4 GeV. 



In addition, the experimental spectra exhibit little de- 
pendence on centrality (except for peripheral collisions), 
and therefore the slope analysis has been performed for 
inclusive qt spectra with dA^ch/dy > 30. In our theoret- 
ical analyses wc focus on semicentral collisions. 

We have tried to follow the same procedure (see also 
Ref. [s^, but as illustrated by the dash-dotted line in 
Fig. [ini a fit to our theoretical spectra (represented by 
the dots) in the above range (0.4 GeV < qr < 1-8 GeV) 
slightly underestimates the effective slopes toward higher 
qx, especially in the lower mass bins. Part of the prob- 
lem is the bias of the fit toward low qx where the yield 
is the largest and thus dominates the total x^- To better 
reproduce the theoretical spectral shape with the above 
function we have therefore chosen to (i) restrict the fit 
range to 1 GeV < qr < 1-8 GeV, and (ii) minimize the 
for the logarithm of the spectra. The resulting fit 
(solid line in Fig. [19]) indeed agrees better with the the- 
oretical spectra. This effect becomes more pronounced 
when the Drell-Yan component is included (which, in 
fact, might indicate that our extrapolation overestimates 
the DY yield at low mass). 

Based on the slightly revised procedure, our slope anal- 
ysis is carried out for the total sum of all calculated qr 
spectra components in twelve AM = 0.1 GeV mass bins 
in the range 0.2 GeV < M < 1.4 GeV (corresponding to 
the experimental ones) for semicentral In-In; by default, 
the Drell-Yan and meson t-channel contributions are not 
included, but studied in separate fits by adding them to 
the totals. As in the experimental analysis, correlated 
open-charm decays arc removed altogether. Our system- 
atic study addresses variations in the EoS as well as the 
possibility of another ~15% increase in the transverse ac- 
celeration of the fireball expansion, i.e., a± = O.lc^/fm 
(which is, in fact, the value we have used in our cal- 
culations for heavy-flavor observables at RHIC [73, as 
extracted from the hydrodynamic model of Ref. [9a|). 
Specifically, we evaluate the following scenarios, as sum- 
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marizcd in the six panels of Fig. 1201 

(a) the default fireball expansion with EoS-A {Tc = 
Tch = 175 MeV), characterized by a thermal freeze- 
out with temperature and radial flow surface veloc- 
ity of (Tfo,/31 fj = (120 McV, 0.57c); 

(b) the default fireball expansion with EoS-B 
(Tc = Teh = 160 MeV) and (Tfo,/31,fJ = 
(136 MeV, 0.57c); 

(c) the default fireball expansion with EoS-C (T^ = 
190 MeV, Tch = 160 MeV) and {Tfo,Pl^fJ = 
(136 MeV, 0.57c); 

(b+) the same as in (b), but with a transverse fireball 
acceleration a± = O.lc^/fm yielding (Tfo,/?^ j^) = 
(135 MeV, 0.65c); 

(c+) the same as in (c), but with a transverse fireball 
acceleration a± = O.lc^/bn yielding (Tfo,/?^ fo) = 
(135 MeV, 0.65c); 

(c]''"'') the same as in (c+), but additionally including the 
contribution from meson-^-channel exchange with 
in-medium p propagator (cf. Sec. IIIDjl : 

(cj^) the same as in (cj^^), but using the vacuum p prop- 
agator in the meson-t-channel exchange contribu- 
tion. 

The upper two panels of Fig. [20] compare the impact 
of the three different EoS using our baseline acceleration 
of aj^ = 0.085c^/fm. The default scenario EoS-A (a) falls 
short of the empirical slopes by a significant margin of up 
to 50-60 MeV in the LMR, while reaching the lower end of 
the data in the IMR. The higher hadronic temperatures 
of scenarios EoS-B (b) and EoS-C (c) improve the situ- 
ation by about 15-20 MeV in the LMR, approximately 
reflecting the increase in the hadronic fireball tempera- 
tures (e.g., Tfo = 136 McV compared to Tfo = 120 MeV 
at thermal freezeout). Scenario EoS-B (b) additionally 
improves around the p peak due to a larger weight of the 
relatively hard components from decays of freezeout and 
primordial p, since the overall thermal hadronic emission 
is smaller than for EoS-A and EoS-C (compare Fig.fTSlto 
Figs. [3 and [HD. However, for both EoS-B and EoS-C the 
slopes in the LMR are still below the data, even though 
the shape is not too bad. 

In an attempt to improve on the slopes, we investi- 
gate a further increase of the transverse acceleration to 
values representative for RHIC energies [tJ HEIj from 
0.085 cVfm to 0.1 c^/fm for EoS-B and EoS-C (scenarios 
(b"*") and (c+), shown in the middle panels of Fig. [^0)1 . 
While this increases the effective slopes at all emission 
stages, it is more efficient for the later hadronic sources 
(in accordance with the blue-shift pocket formulae given 
in Sec. IVBp . Consequently, the LMR exhibits further 
improvement, and, importantly, the consistency with the 
data in the IMR is not spoiled due to the fact that the 



main contributions there arise from earlier stages with 
higher temperatures and less flow (dominated either by 
the contributions from the QGP phase for EoS-B or four- 
pion annihilation in the high-density hadronic phase for 
EoS-C, cf. Fig. [18]). At this point, EoS-B provides some- 
what better agreement and one might be tempted to 
take this as evidence for a prevalently partonic emission 
source, as suggested in Refs. |87ll88|. However, remaining 
uncertainties preclude an unambiguous interpretation, as 
we will now show. 

In the two bottom panels of Fig. [20] we focus on the 
EoS-C scenario with a± = O.lc^/fm. The left panel illus- 
trates the sensitivity of the slopes to the additional inclu- 
sion of i-channel meson exchange processes. As discussed 
in the previous section, their spectra alone carry a slope 
very similar to the data; not surprisingly, their inclusion 
increases the slopes in the LMR with no impact on the 
IMR slopes. When implementing the t-channel meson 
exchanges with the vacuum p propagator, the effect on 
the slopes is somewhat larger compared to using the in- 
medium Dp, especially in the free p mass region. A sin- 
gle np scattering with unmodified p mesons would most 
closely represent interactions of p mesons in a kinetic 
regime (rescattering without thcrmalization) , and thus 
provide a missing link between the two extremes of ther- 
mal radiation and primordial p's surviving jet quenching 
(no reinteraction) in our baseline calculations (note that 
there is no contribution from "quenched" p's). The over- 
all level of agreement with the data is now very similar 
to EoS-B. Also recah that Tc = 190 McV maximizes the 
ratio of hadronic to QGP emission; for Tc = 175 MeV 
four-pion emission still outshines the QGP in the IMR, 
while the slope parameters are closer to the EoS-B sce- 
nario. The effect of the Drell-Yan contribution on the 
slopes is shown in the bottom-right panel of Fig. [^D] In 
the IMR the change in slope is consistent with the experi- 
mental results (full vs. open squares), while in the lowest 
mass bins the uncertainties in the extrapolation of the 
DY are large. We emphasize again that more stringent 
constraints can be enforced once an accurate knowledge 
of the absolute normalization in the experimental mass 
and momentum spectra becomes available. 

There are uncertainties in the experimental extraction 
of the slope parameters beyond the statistical error bars 
shown in the figures. In Ref. (ssj a systematic error of 
~7 MeV is estimated due to the subtraction of back- 
ground and decay (cocktail) sources. For the LMR, it 
has been estimated that the subtraction of a potential 
Drell-Yan contribution could lower the effective slopes 
by 5-10 MeV. One should also be aware of that the ap- 
plicability of hydrodynamic (or thermal fireball) descrip- 
tions is limited in transverse momentum, typically to 
Pt < 2 GeV for pion spectra at the SPS (consistent with 
where the spectra reach the collision-scaling regime in 
Fig. [6]), and to significantly smaller px for elliptic flow. 
Fireball models might compensate for this to some ex- 
tent by implementing a larger acceleration, which may 
be (part of) the reason that the dilcpton slopes prefer 
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FIG. 20: (Color online) Effective sfope parameters fitted to the theoretical qT spectra using Eq. (|45|l for the following scenarios 
(by default, DY and f-channel meson exchange reactions are not included): top left: EoS-A (dashed line) vs. EoS-B (dash-dotted 
line), both with ax = 0.085c'^/fm; top right: EoS-A (dashed line) vs. EoS-C (dash-dotted line), both with ax ~ 0.085c'^/fm; 
middle left: EoS-B with either ax = 0.085c^/fm (dash-dotted line) or ax ~ O.lc^/fm (dashed line); middle right: EoS-C with 
either ax = 0.085c^/fm (dash-dotted line) or ax = O.lc^/fm (dashed line); bottom left: EOS-C with ax ~ O.lc^/fm (lower 
dashed line), additionally including meson f-channel contributions with free (dash-dotted line) or in- medium (middle dashed 
line) p propagator; bottom right: EOS-C with ax = O.lc^/fm (dashed line), additionally including Drell-Yan annihilation 
(dash-dotted line). 



aj^ = O.lc^/fm. To further scrutinize the spectra, 
it would be very illuminating to reiterate earlier imple- 
mentations of the dilepton rates from hadronic many- 
body theory into hydrodynamic [sol (w ith a more real- 
istic hadro-chemsitry) or transport [49| simulations. In 



this context, a measurement of dilepton elliptic flow could 
provide further insights in the detailed emission history 
of the medium at SPS energies [90| . 
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E. Update of Comparison to CERES/NA45 Data 

After the refinements in our theoretical approach for 
SPS dilepton spectra in the context of the NA60 data, 
we revisit in this section the consequences for the com- 
parison with CERES/NA45 in semicentral and central 
Pb-Au collisions, which showed good agreement in our 
earlier wor ks [3, m, [e!]. In the following, we con- 
strain ourselves to the above default scenario for the fire- 
ball evolution with a transverse acceleration of aj_ = 
0.085 c^/fm. As in the In-In case, the increased ex- 
pansion rate entails shorter fireball lifetimes; we choose 
the latter as to obtain the same thermal freezeout tem- 
perature for central (Tfo ~ 106 MeV) and semicentral 
(Tfo ~ 112 MeV) Pb-Au collisions as in our previous 
work [IJl, implying a reduction by 2 fm/c for both 
centralities (13.5 11.5 fm/c and 12 10 fm/c, re- 
spectively). Also, since for the larger system sizes {i.e., 
lifetimes) the freezeout- and cocktail-p contributions [1| 
are significantly suppressed relative to thermal radiation, 
and since we are not interested in (high-) qt spectra here, 
we adopt the (old) simplified treatment for the freeze- 
out p b y r unning the fireball an additional duration of 
1 fm/c [2J|. The larger lifetimes relative to In-In colli- 
sions furthermore imply that the nonrelativistic formula 
for the surface expansion velocity, = a±t, approaches 
values uncomfortably close to c in the late stages of the 
collision. To implement a relativistically covariant accel- 
eration in a simple way, we model the acceleration of the 
fireball in analogy to the relativistic motion for a charged 
particle in a homogeneous electric field. The radial sur- 
face velocity and fireball radius then take the form 



a^t/c 



^2 



r = ro + 



(46) 



which tames the acceleration at late times and matches 
the non-relativistic case for early times. 

In the upper panel of Fig. [21] the updated calculations 
for dielectron spectra in semicentral Pb-Au collisions are 
compared to our previous results (using the same hadro- 
chemistry and p spectral function) , including experimen- 
tal acceptance cuts as defined by CERES/NA45. For 
masses M > 0.6 GeV, the p spectral function yield is re- 
duced by 15-20%, directly reflecting the reduction in the 
lifetime of the hadronic phase in the fireball. At lower 
masses, however, the difference is smaller and eventually 
disappears for M < 0.35 GeV. The reason for this is the 
interplay of the single-electron cuts and the increase in 
transverse flow, where the latter gives larger momentum 
to both electrons which in turn have a larger probability 
to make it into the acceptance. For the QGP contribu- 
tion, the difference is smaller (less than 10%) and weakly 
dependent on mass, since the increased transverse flow 
plays little role in the early phases (we recall that, for 
the canonical choice of the formation time, tq = 1 fm/c. 
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FIG. 21: (Color online) Upper panel; comparison of ther- 
mal dilepton spectra in semicentral Pb-Au collisions (includ- 
ing experimental acceptance) for p (short-dashed lines) and 
QGP emission (dash-dotted lines) using our previous (labeled 
"old") and updated fireball evolution. The solid line is the 
sum of in-medium p, ui and ((> contributions, while the long- 
dashed line results from four-pion annihilation including chi- 
ral mixing as described in Sec. Ill CI The Drell-Yan contri- 
bution is below the displayed range; correlated charm decays 
are neglected. Lower panel: The solid line is the total in- 
medium thermal emission -I- hadronic cocktail (dash-dotted 
line, adopted from the experimental evaluation H]), compared 
to combined 1995/1996 CERES/NA45 data [JT The short- 
dashed line is the sum of thermal emission with a vacuum p 
spectral function and the cocktail. 



the initial temperature for 30% central Pb-Au amounts 
to To = 203 MeV, and 210 MeV for 7% central). The 
addition of in-medium uj and decays affects the to- 
tal yield mostly around the respective free masses, more 
significantly in the region. The in-medium four-pion 
contribution becomes relevant for masses AI > 1 GeV, 
with a relative strength compared to the other sources 
that closely resembles the In-In case. The Drell-Yan con- 
tribution docs not show up on the scale of Fig. [^T] (and 
has been restricted to masses above '--^l GeV). Correlated 
charm decays are not accounted for. The lower panel of 
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FIG. 22: (Color online) Comparison of theoretical calcula- 
tions of dilepton excess spectra in central Pb-Au collision, 
compared to CERES /NA45 data [H] ; for the line labeled "to- 
tal (no bar p)" the full in-medium p spectral function has been 
replaced by the one which only includes medium effects in a 
meson gas. 

Fig. [2T] summarizes the comparison of the theoretical cal- 
culations with the combined 1996/1996 CERES/NA45 
data including the hadron decay cocktail. The level 
of agreement is very similar to the one with our earlier 
results 11,163. 

Finally, we turn to the most recent CERES data 
for central Pb(158 AGeV)-Au collisions ^], which are 
shown in Fig. [51] in a form similar to the NA60 data in 
the previous sections, i.e., on a linear scale and with the 
cocktail subtracted. Again, our calculations employing 
the updated fireball describe these data fairly well. Vari- 
ations of ±20% in the fireball lifetime (affecting the total 
dilepton yields at the same level with little change in the 
spectral shape), would still lead to reasonable agreement. 
The importance of baryon- induced medium modifications 
at SPS energies is reiterated by the dash-double-dotted 
curve where only meson-gas effects are included in the 
p spectral function; this scenario does not properly re- 
produce the experimental spectra. When comparing to 
the NA60 dimuon mass spectra, the four-pion and QGP 
contributions in the IMR appear smaller relative to the 
signal in the mass region below the free p mass. This is 
essentially due to the experimental acceptance, which for 
CERES/NA45 is significantly larger in the low-mass re- 
gion than for NA60, especially for low-momentum lepton 
pairs. 

A particularly striking feature of the theoretical predic- 
tions is the large enhancement below the two-pion thresh- 
old, which is closely connected with the approach to (and 
constraints from) the photon point. The very-low-mass 
excess is rather sensitive to the baryon-induced medium 
effects as well, as illustrated by the curve where the lat- 
ter are switched off. Experimentally, this regime is only 
accessible with dielectrons (and after subtraction of the 



7r° Dalitz decay), and the first hint from data is consis- 
tent with the theoretical yield. Clearly, it would be very 
worthwhile to further explore this region. 



VI. SUMMARY AND CONCLUSIONS 

In the present paper, we have conducted a quantitative 
study of dilepton spectra in semi-/central heavy- ion colli- 
sions at full SPS energy {^/s = 17.3 GeV). We have sup- 
plemented our earlier calculations of thermal radiation 
(which reproduce available mass spectra) with sources 
of non-thermal origin (expected to become relevant at 
high transverse momentum and for small systems), and 
we have scrutinized the results with respect to hadro- 
chemical and flow properties of the underlying thermal 
fireball expansion. 

The key quantity to describe thermal dilepton rates 
is the (imaginary part of the) electromagnetic current 
correlator. In hadronic matter, we have evaluated its 
medium modifications in terms of many-body spectral 
functions for the light vector mesons at low mass, and chi- 
ral vector-axialvector mixing for the continuum part at 
intermediate mass. Partonic emission above Tc has been 
approximated by perturbative quark-antiquark annihila- 
tion. An appealing feature of this description is a smooth 
merging of the hadronic rates with the partonic ones 
at temperatures around the expected phase transition, 
Tc ~ 160-190 MeV. As a new ingredient we have com- 
puted thermal dilepton rates induced by <-channel meson 
exchange in 7rp — > tt/"*"/" reactions, which we found to 
become relevant at momenta above ~1.5 GeV. We have 
augmented our approach by non-thermal sources, i.e., 
Drell-Yan annihilation, as well as decays of p mesons 
from primordial production (subject to jet quenching) 
and from decoupling at thermal freezeout. We have em- 
phasized that these p decays carry an extra kinematic 
Lorentz-7 factor relative to thermal radiation, resulting 
in harder dilepton-^/T spectra. 

All sources have been implemented into a space-time 
evolution for A-A collisions at SPS, for which we em- 
ployed an isentropic fireball model with QGP, mixed and 
hadron-gas phases and updated expansion parameters to 
better match the absolute dilepton yields measured by 
NA60. Hadrochemical freezeout and subsequent evolu- 
tion are constructed in line with measured hadron multi- 
plicities. The calculated dimuon invariant-mass spectra 
agree well with NA60 data in central and semicentral In- 
In collisions, including qT bins below 0.5 GeV and above 
1.0 GeV, thus confirming our earlier results for thermal 
emission and validating the predicted in-medium vector 
spectral functions (dominated, but not exhausted, by the 
p). When comparing to qx spectra, the new sources 
contribute to provide good agreement for central colli- 
sions, while discrepancies persist with semicentral data at 
(7t > 1 GeV for masses around and below the free p mass. 
In particular, effective slope parameters for gr > 1 GeV 
amount to T^s ~ 160-210 MeV, --25% short of the em- 
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pirically extracted values in the LMR. 

We have investigated the sensitivity of our results to 
variations in the critical and chemical-freezeout temper- 
ature of the fireball. The spectral shape of the invariant- 
mass spectra turned out to be insensitive for the range 
Tc = 160 - 190 MeV. The reason is the "duality" of 
the thermal emission rates around Tc, rendering the 
"melting" of the p meson a robust signal. A low value 
of Tc ~ 160 MeV entails that QGP radiation domi- 
nates over hadronic emission in the IMR. A low value of 
Teh = 160 MeV implies smaller pion chemical potentials 
(and thus higher temperatures) in the hadronic evolu- 
tion, increasing the effective slope parameters by about 
15-20 MeV, still -30 McV short of the data in the LMR. 

In an attempt to resolve this discrepancy, we have 
implemented a ~15% increase in the transverse fireball 
acceleration as previously employed for central Au-Au 
collisions at RHIC energy. This elevates the slope to 
^^250 MeV around the free p mass (and properly re- 
produces the increase below) , while still being consistent 
with the (upper) experimental value of 205 ± 20 MeV in 
the IMR where the blue shift effect is less pronounced due 
to the prevalently early emission. In particular, we have 
demonstrated that, within the current theoretical and 
experimental uncertainties, both scenarios for a predom- 
inant emission source in the IMR - hadronic (four-pion) 
or partonic - are viable. In either case, the radiation 
mostly emanates from matter at temperatures around 
Tc ~ 160 — 190 MeV. An unambiguous distinction be- 
tween a partonic and hadronic source in the IMR there- 
fore appears difficult at SPS energies. The larger trans- 
verse flow developed in the partonic stage at RHIC en- 
ergies could facilitate this task [9l|, although additional 
complications arise due to the much larger open-charm 
contribution and its in- medium modifications. 

Finally we have checked that our refined assessment of 
dilepton sources, together with the improvements in the 
fireball expansion, preserves our previously found agree- 
ment with the CERES/NA45 diclcctron spectra in Pb-Au 
collisions. The larger system size leads to an apprecia- 
ble increase in thermal radiation, thus reducing the un- 



certainties associated with the modeling of the break-up 
stage. In addition, dielectrons allow access to the mass 
region below the two-pion threshold, where the medium 
effects on the p-spectral function, augmented by the ther- 
mal Bosc factor (and photon propagator), predict a large 
thermal dilepton signal. First data in this region support 
this signal, and further corroborate the importance of 
baryonic interactions for the in-medium p spectral func- 
tion. 

Future efforts should be pursued along several direc- 
tions: first and foremost, the theoretically calculated 
vector spectral functions must be extended to their chi- 
ral partners (most notably the axialvector (ai) channel 
as the partner of the p), to establish direct connections 
to chiral symmetry restoration. The evaluation of in- 
medium Weinberg sum rules (9^ , in connection with con- 
straints from lattice QCD, will play an important role in 
this enterprise [9l, [l^l- Second, from the phenomenolog- 
ical side, the implementation of emission rates into both 
hydrodynamic and transport simulations should be re- 
iterated, to compute qt spectra and check whether the 
rather large acceleration of the fireball model can be justi- 
fied for the In-In system at SPS. The study of dilepton el- 
liptic fiow [sOl could shed further light on the time profile 
of the radiation. At RHIC energies, partonic collectivity 
increases significantly; an unexpectedly large dilepton ex- 
cess in the LMR recently observed by PHENIX needs to 
be understood [9^ . Work in several of these directions is 
in progress. 
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